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ABSTRACT 

In this thesis, the calculation of the electric dipole 
polarizability tensor of the He 2 dimer is described, and the results 
are used in the computation of several dielectric and optical 
properties of helium gas, at both high (322 0 K) and low (4°K) 
temperatures. The properties considered are the second dielectric 
virial coefficient, the second Kerr virial coefficient, and the 
depolarization ratio of the integrated intensities for the Raman 
scattering experiments. 

The thesis consists of five parts. In the first part, the 
polarizability and various properties are defined. In the second 
part, the calculation of the polarizability in the long-range region 
in terms of a quantum mechanical multipole expansion is described. 
The formulas which are obtained are applied to the H 2 , He 2 , and HeH 
diatomics, and the results are compared to those of model 
calculations. In the third part, the calculation of the He 2 
* 
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polarizability in the overlap region via coupled Hartree-Fock 
perturbation theory is described, and a basis set selection 
procedure is delineated which allows the Hartree-Fock limit to be 
approached with sets of reasonable size. It is further shown that 
the long-range limit of the coupled Hartree-Fock polarizability is 
identical to that of the point dipole model. The calculation of the 
quantum pair distribution function for both the 3 He and 4 He isotopes 
at 4°K is discussed in the fourth part. The calculated values of 
the properties of helium gas are given in the fifth part. The order 
of the computed second dielectric virial coefficients for 3 He and 
4 He at 4°K is found to be opposite to that of experiment. This 
experimental -computational discrepancy is discussed in detail; while 
there is an uncertainty in the computed results due to the neglect 
of electron correlation, it is concluded that the discrepancy is 
probably due to inaccuracies in the experiments. 
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CHAPTER I 
INTRODUCTION 


A. Motivation 

"We are perhaps not far removed from the time when we 
shall be able to submit the bulk of chemical phenomena 
to calculation." (J. L. Gay-Lussac, Memoires de la 
Societe' d'Archeil, 2, 207 (1808).) 

"It is also a good rule not to put overmuch confidence 
in the observational results that are put forward until 
they have been confirmed by theory." (A. Eddington in 
The Coming of the Golden Aoe by G. S. Stent.) 

There is a well established paradigm for computing the macroscopic 

properties of the rare gases. The first step in the calculation involves 

treating the electronic structure of the constituent atoms or molecules 

of the gas via quantum mechanics. The link between the molecular 

properties and experiment is provided by statistical mechanics. The need 

for either classical or quantum statistical procedures is dictated by the 

parameters of the particular experiment which one is attempting to 

describe; low temperature experiments on light atoms, for example, 

necessitate the inclusion of quantum statistical effects. Then, since 

statistical formulas which relate experimentally observable quantities to 

the properties of the molecular system are available, one can, at least 

in principle, obtain results for the macroscopic system using a 

completely computational procedure. 
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In this light, experiments are an alternative to computation. 
While the approach actually taken in a given case is largely determined 
by convenience, the computational scheme can result in a more detailed 
understanding of the nature of the physical interactions which give 
rise to the phenomenon under investigation. 

In the work which is reported in this thesis, we adopt the 
computational approach, and calculate a number of dielectric and 
optical properties of gaseous helium. 

The quantities in which we are interested are: the pressure 

dependence of the dielectric constant; 1 the static electric-field- 
induced birefringence (Kerr effect); 2 and the collision induced Raman 
spectrum. 3 The experimental results, as well as the theory, are 
formulated in terms of virial expansions in powers of the density. We 
have considered the terms which are linear in the density and arise 
from interactions involving pairs of helium atoms. The calculation 
of each of these properties involves both the electric polarizability 
tensor (quantum mechanics) and the quantum radial distribution function 
for helium (statistical mechanics). The corresponding properties for 
various systems in addition to helium, including other rare gases, 
have been experimentally studied by other workers. We, however, 
restricted our work to helium because it is the only system which 

currently may be treated accurately using a completely non-empirical 
approach. 
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There are several reasons for undertaking this study, not the 
least of which is the fact that the purely computational procedure has 
never before been carried to completion (for helium) for a property 
which requires the accurate calculation of both a molecular property 
(other than the energy) and the quantum distribution function. Even 
for the relatively simple case of helium, approximations are made in 
both the quantum mechanical and statistical mechanical parts of the 
computation. While it is not our intent to question the general 
approach, we wish to examine the sensitivity of the experimental 
properties to these approximations, and to assess the interpretation 
of the experimental data in terms of molecular properties. One hopes 
that a critical comparison of experimental and computed results would 
either verify the utility of a given set of approximations or indicate 
the direction in which one should move in order to improve the 
approximations . 

Most of our work concentrates on calculating how the polarizabil ity 
of a pair of helium atoms changes as a function of internuclear 
separation. For collisions at thermal energies, these changes probably 
amount to no more than a few percent of the sum of the polarizabilities 
of the isolated helium atoms. Nevertheless, these changes give rise to 
the observable effects mentioned above. The accurate calculation of 
these changes is a difficult computational problem. 

Quantum mechanical calculations of molecular polarizabilities 
have had a long history, dating from Schrodinger's second paper wherein 
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the polarizability of the hydrogen atom was obtained . 4 For more 
complicated systems than this, however, polarizability calculations 
are far from routine. Although semi -empirical methods have been 
developed which are capable of giving good upper and lower bounds to 
the exact polarizability , 5 most non-empirical calculations are for 
atoms, or molecules at fixed geometry, and are at the level of the 
independent particle {Hartree-Fock) model. Only recently have atomic 
polarizability calculations, for systems other than helium, been 
reported in which electron correlation is taken into account . 6 For 
atomic helium itself, the calculations of Buckingham and Hibbard have 
yielded the polarizability to six significant figures . 7 For molecules, 
only the fully correlated calculations of Kolos and Wolniewicz on H 2 
can be considered definitive . 8 

To anticipate the pattern of our calculation of the helium pair 
polarizability, it is useful to briefly recall the history of 
calculations of the pair potential . It has been known since Lennard- 
Jones' successful fit of equation of state and transport data 8 that 
the helium pair potential is characterized by long range attraction 
and short range repulsion, resulting in a shallow van der Waals well 
at intermediate separations. Until recently, there were two separate 
computational approaches to explain this behavior: the multipole 
expansion of the potential at large internuclear separations and the 
Hartree-Fock approximation at shorter separations. The two approaches 
make different assumptions: the former assumes that the overlap 
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between the atomic wavefunctions is negligible and treats the 
interatomic correlation which gives rise to the attractive multipole 
interactions; the latter neglects correlation but allows the atoms 
to overlap, resulting in a purely repulsive interaction. It has 
only been within the last several years that calculations have been 
performed which employ a wavefunction which passes smoothly from 
long range to the region of the van der Waals well. 10 Calculations 
of the helium pair potential are now in good accord with experimental 
data, although there are still some questions regarding the true depth 
of the well . 10b * 11 

In our calculations of the pair polarizability we have followed 
the pattern of the pair potential calculations, and have considered 
both the long range multipole expansion and the coupled Hartree-Fock 
approximation. As we shall show, in contrast to the pair potential 
case, the coupled Hartree-Fock model for the pair polarizability does 
reduce to the multipole expansion at large separations, although the 
multipole coefficients are incorrect due to the neglect of electron 
correlation. Thus, there is an uncertainty in our final results for 
the pair polarizability, which can only be removed by using a 
wavefunction which accurately describes electron correlations. 

Our statistical mechanical calculation of the radial distribution 
function is more reliable. The essential ingredient here is an accurate 
pair potential function. We have used the MDD-2 potential 12 which fits 
well the available experimental data on macroscopic and scattering 
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properties of helium and also agrees with the non-empirical 
calculations. Further, we shall show that the dielectric and optical 
properties are not as sensitive to the pair potential as they are to 
the pair polarizability. In the calculation of the quantum 
distribution function, we followed methods established recently by 
Klemm and Storer. 13 

In the next sections, we will review the experimental and prior 
theoretical work on the pair polarizability of helium and on the 
aforementioned dielectric and optical properties, and summarize our 
own computations. Later chapters contain detailed explanations of 
the different aspects of our work. 

B. Polarizability Definitions 

The polarizability of an atom or molecule is a measure of how 
easily its electronic charge cloud is distorted by an external electric 
field. The electronic structure in the absence of any external forces 
is determined primarily by the balance of the electron-nuclear 
attractive forces with the electron-electron and nuclear-nuclear 
repulsive forces, in accord with the Schrodinger equation and the 
symmetry requirements of the Pauli Principle. The measure of these 
internal atomic fields is the ratio of the electronic charge to the 
square of the Bohr radius, e/a^ - 10 9 volts/cm . Since most 
laboratory fields are considerably less than 10 5 volts/cm , the 
effect of such applied fields on the electronic structure of the 
molecule is expected to be small. 
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For applied electric fields which are sufficiently small, it is 
natural to extend the field-free treatment of molecular electronic 
structure to the field-dependent case via perturbation methods. 

Thus, one assumes that the electronic energy, E , of a molecule in 
an external electric field, F = (F V ,F ,F'J can be expanded in 
powers of the components of F : 


E = E<°> - I „j°> F, - M Z. a 1 jFl F. 'HI 


1 >J 


i > j . k 


( 1 . 1 ) 


" Y WA F 1 F j F k F * + '** 


In this equation, is the i-th component of the permanent 

dipole moment of the molecule, a., is the ij component of the 

* J 

polarizability tensor, while the hyperpolarizabilities 6 and y 
refer to tensors of rank greater than 2 and reflect higher field 
effects. The external field also affects the total dipole moment 
of the molecule, and the i-th component of the dipole moment is 
given by 


_ / 3 E \ 

y i l 3F. ; F.=0 ’ 


( 1 . 2 ) 


Then, from Eq. (1.1) for the energy, 


"i s u ! 0> + 1 “i/j + 


(1.3) 
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These last two equations suggest two ways by which the polarizability 
can be obtained. That is, it can be obtained from a second order 
expression which involves the energy, 


U 


/2 'sF,3F, ; F=0 > 

■ J 


0.4) 


or, from a first order expression which involves the dipole moment. 


3p.j 

a ij = ^F7^F=0 • 

J ' 


(1.5) 


From Eq. (1.1) it is seen that the effect of nonzero 
hyperpolarizability terms is to give rise to deviations from a linear 
polarization law. Buckingham has written an excellent review article 
on induced electric moments and their modification by intermolecular 
interactions. 14 

To place the above definitions in a quantum mechanical framework, 
we consider the Schrodinger equation for a molecule in a constant 
electric field: 


[H ( F) - E(F)MF) = 0 . (1.6) 

In this equation, H(F) is the field-dependent Hamiltonian operator 
for a single molecule, and 'i'(F) and E(F) are the (field-dependent) 
wavefunction and energy respectively, for the molecule. The exact 
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solution of this equation is generally impossible to obtain, so that 
one resorts to approximate schemes. For molecules, one begins by 
making the Born-Oppenheimer approximation, and then applies either 
perturbation theory or variational methods to obtain approximate 
solutions to the equation which results. 

For constant field, the Born-Oppenheimer Hamiltonian for a 
molecule having N electrons whose positions are denoted by r. , 
and M nuclei of charge Z a at R a , can be written as: 



(1.7) 


Here, u is the molecular dipole moment operator defined (in atomic 
units) by: 


A 

u 


N 

■ l :i 
i-1 1 


M 

I Z R , 

a-1 a ~ a 


0 . 8 ) 


and H 0 is the field-free, Born-Oppenheimer Hamiltonian for the 
molecule. The essence of the perturbation method is that one assumes 
that the ground state energy, Eg(F) , and the corresponding 
wavefunction , ^(F) , can be expanded in a power series in F = F f 


E 0 (F) = l F n E<") 
n=0 


» 


(1.9) 
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and 

¥o(F) = l F n . (1.10) 

n-0 

Substituting these last two equations into Eq. (1.6), and setting 
the coefficient of each power of F equal to zero, leads to the 
usual set of perturbation equations. The zero and first order 
equations are respectively: 

(H 0 - - 0 , (1.11) 

and 
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where we have used the normalization condition 


+ <*o 0) l v o 1}> = 0 


(1.14) 


Comparing Eqs. (1.1) and (1.9) we see that the quantum mechanical 
analogue of Eq. (1 .4) is 

\ 3 = - JC 0 2) F2 • 0 - 15 ) 

or. 

I - 2<4° ) i5-f|f< 1) >F2 . (1.16) 

To obtain the analogue of Eq. (1.5), we note that quantum 
mechanically the dipole moment is obtained as an expectation value: 

y * <u> = <¥(F)|jj|v(F)» . (1.17) 

In view of Eqs. (1.9) and (1.10), this expectation value can also be 
represented as a power series in the field: 


y 


l l F 

n»0 m=0 


n+m 




(1.18) 


= <T 


( 0 ) r . ( 0 ) 
n 1^ Kn 


> + 2<'4° ) !y|^ 1 )>F + .... 


(1.19) 
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Comparing this last equation to Eq. (1.3) we see, of course* that 
the same definition of a (Eq. (1.16)) results from consideration 
of either the energy or the dipole moment within the context of 
perturbation theory. 

The problem of computing the polarizability has been reduced to 
that of solving Eq. (1.12) for . Using a reduced resolvent 

notation, R 0 = (E^ - H 0 ) _1 , one writes the solution of Eq. (1.12) 
formally as : 

= “R 0 ?-f Y o 0) * (1.20) 

so that a , the polarizability tensor, becomes : 

a = -2 <t( 0) |u R 0 jj|*J 0 )> . (1.21) 

This last equation is the immediate starting point for the 
multipole expansion of a which will be considered in the next chapter. 

If one uses variational methods rather than perturbation theory 
to solve Eq. (1.6), then Eq. (1.21) is not the expression for the 
polarizability which is used. This is our situation in the overlap 
region. Here we optimize the variational functional 

E = <4'|H(F)|y>/<y|^> , (1.22) 
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to obtain a field-dependent; trial wavefunction ¥ ( F) , from which 
the dipole moment is approximated as 

= <4'(F)[p i |y(F)>/<*|¥> (1.23) 


To calculate the polarizability we return to Eq. (1,5) and use the 
definition of derivative explicitly: 


a . 




1 ini 

F -► 0 



(1.24) 


This is the starting point of the Rartree-Fock calculations which 
will be discussed in Chapter HI. 

With respect to calculations in the short range region, one 
could, from a strictly formal point of view, just as easily use the 
definition of the polarizability in terms of the energy, Eq. (1.4). 
However, from a computational standpoint, the precision in the 
difference between the field-dependent and field-free energies is 
less than the corresponding dipole moment difference because one 
uses small applied fields (to eliminate hyperpolarizability terms) 
and the polarizability is a second order effect in the energy while 
it is a first order effect in the dipole moment. 

The experimental quantities which will be given in the next 
section involve a knowledge of the polarizability tensor for at most 
a pair of interacting helium atoms. If, for a cylindrical ly 
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symmetric system such as this, we choose the z-axis to lie along the 
internuclear axis, the polarizability tensor is diagonal, and Eq. 
(1.24) becomes: 


r !I(_L) 


lim 


y iiu) 


- u (0) 

11 ( 1 ) 


HD 


'II(-L) 


(1.25) 


where the subscripts B and ± refer respectively to directions 
along and perpendicular to the internuclear axis. For the case of 
interest In this work, the term involving is absent because 

the field-free dipole moment of a homopolar diatom is zero. 


In terms of the components given in Eq. (1.25), we define the 
mean polarizability of the diatom (internuclear separation R ) and 
the anisotropy respectively as follows: 


a(R) = (a,(R) + ?«j_(R))/3 

0.26) 

B(R) = a ||(R) - cij^R) . 

(1.27) 


We further define the mean incremental polarizability by: 

Ace(R) = ot(R) - 2a 0 , (1.28) 

where ct 0 is the polarizability of an isolated helium atom. Note 
that because of the symmetry, a knowledge of the quantities a and 



g is equivalent to a knowledge of aj and aj_ . Hence, the 
measureable quantities discussed in the next sections provide 
information about the individual tensor elements. 

We now turn to a consideration of the experiments which measure 
the polarizabilities Act and b . 

C. Pressure Dependence of the Dielectric Constant 

When a rare gas sample is subjected to an external electric 
field, the value of the total field at a point in the sample is not 
equal to the applied field, because of additional fields arising 
from the (induced) dipole moments of the atoms. In the classical 
treatment of this problem, Lorentz in 1908 approximated the local 
field in terms of molecular polarizabilities and obtained the 
Clausius-Mosotti formula for the dielectric constant, 15 

V m = tt N a 0 . (1.29) 

Here, V m is the molar volume, N is Avogadro's number, e is the 

dielectric constant, and a 0 is the polarizability of an isolated 

atom. The left-hand side of this equation is known as the Clausius- 

Mosotti function. 

Experimental determination of the Clausius-Mosotti function 
requires the measurement of the dielectric constant of the gas at a 
given temperature and pressure. This is done by measuring the ratio 
of the capacitance of a sample cell to that of an evacuated cell. 
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Early experiments by Uhlig, Kirkwood, and Keyes showed that the 
Clausius-Mosotti function exhibits a dependence on the gas density. 16 
In an attempt to explain the discrepancy between Eq. (1.29) and 
experiment, Kirkwood noted that the following two assumptions which 
are made in the derivation of Eq. (1.2S) are invalid: 17 

1. The polarizability is itself independent of density; 

2. The dipole moment of an individual ato?ns does not change 
as the atom-pairs move through their various phases of 
thermal motion. 

In reformulating the problem without these assumptions, Jansen 
and Mazur 10 and Buckingham and Pople 19 made use of classical 
statistical mechanics to relate the dielectric constant to the total 
induced dipole moment of the system, 


where the 


e-1 

e+2 


m 


electric moment 



<M(F)> 


n. dwc 

F -> 0 Y 
is given by: 


(1.30) 


<M(F)> 


•f C!3 (t>F) ♦fjexp[-(U-M(x ,F) »F)/kT]d T 
/exp[-(U-M(i,F) •F)/kT]dx 


(1.31) 


Here, M(x,F) is the total electric moment of the system, when the 

molecular configuration is T , and the applied field is F . It is 

N 

assumed that M = £ y where y. is the dipole moment of the k-th 

k=l ~ K 



17 


molecule, and N is the number of molecules in the specimen. In 
addition, in Eq. (1.31), f is a unit vector in the direction of F , 
and U. is the intermolecuiar potential. 

By making a virial expansion of Eq, (1.30), and applying the 
results to S-state atoms, which have neither a permanent nor a 
collision-induced dipole moment, they obtained: 


e-1 

e+2 


m 


A + B V" 1 + C V- 2 + 
e e m e m 


(1.32) 


where 


and. 


A 


£ 


4uN 

— a G , 


(1.33) 


B e = iliNif dRR 2 MR)e -U(R)/kT _ 
0 


0.34) 


In^Eq. (1.34) the Aci(R) is the incremental mean polarizability of 
a pair of atoms with separation R defined in Eq. (1.28), k is 
Boltzmann's constant, and T is the absolute temperature. The 
coefficients A and B are called dielectric virial coefficients. 
McQuarrie and Levine extended these results to include the derivation 
of, an expression for the third virial coefficient, C , for the 

G 

case of non-polar molecules with axial symmetry, such as H 2 or C0 2 . 20 
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The A e and B £ were measured for several rare gases at room 
temperature by Orcutt and Cole. la In addition, low temperature 
data (T = 4°K) for both 3 He and ‘♦He were taken by Kerr ana 
Sherman. lb We present the results of these experiments in Table 1.1. 

The percent change in the Clausius-Mosotti function at room 
temperature is on the order of .1% for helium and 1% for argcn for 
pressures near 2 atm. At low temperatures (near 4°K) , the change 
is some 20 % less than the room temperature change. Thus for 
pressures on the order of 1 atm. or less, it is very difficult to 
measure the pressure effects on the Clausius-Mosotti function. The 
helium date at 4°K of Kerr and Sherman exhibit a large amount of 
scatter. This is because the inherently small effect is compounded 
by the small amount of sample present in the cell at these pressures. 
The amount of scatter is reduced at higher pressures (~4 atm.}, but 
then the validity of the virial expansion is doubtful. 

Several workers have coupled the results of model polarizability 
calculations with classical statistical mechanics to obtain high 
temperature estimates of B £ . ib. 19 *2i Some of the earliest of these 
calculations employed the point dipole approximation. 19 The 
simplicity of this approximation makes it particularly attractive, 
hence we now treat it in more detail. With this model, one assumes 
that the dipole moment induced in a given atom (denoted by 1) due to 
action of the external field F , and a second atom (denoted by 2), 
a distance R away, is given by: 
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Table 1.1. Experimental first and second dielectric virial 

a 

coefficients for several rare gases. 


System 

T(°K) 

A £ (cm 3 /mole) 

• ■■ .■ • ■=.=■•= — - ■ 

B (cm 6 /mole 2 ) 
e 

He 

322 

0.519 

± 0.001 

-0.06 

± 

0.04 

Ne 

322 

0.998 

± 0.001 

-0,03 

± 

0.10 

Ar 

322 

4.142 

± 0.002 

0.39 

± 

0.20 

Kr 

322 

6.267 

± 0.003 

5.6 

+ 

0.30 

3 He 

4 

0.516951 

± 0.000088 

-0.030 

± 

0.004 

4 He 

4 


— 

-0.023 

± 

? 


a The room temperature results are those of Orcutt and Cole, 
reference (la); the low temperature results are those of Kerr and 
Sherman, reference (lb). 
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vji = a o(f + Fi) , 


(1.35) 


where F L is the field at 1 due to the dipole at 2 . A similar 
expression is used for y 2 . Using the usual expressions for the 
field due to a dipole, it is seen that yj can be written: 


F x a o F y a o F z 
M+a' 0 R-3 * l+a 0 R“3 * l-2a 0 R-3 ; 


(1.36) 


where a 0 is the polarizability of an isolated atom. Differentiating 
Eq.' (1.36) with respect to the components of F » and expanding the 
result for large R , yields: 


a i (R) = 2a 0 + 4aJ R" 3 + 8aJ R“6 + _ j (1.37) 


and 


a ± (R) = 2a 0 - 2a* R-3 + 2 «q R" 6 + .... (1.38) 

These results are exact through order 3 in R" 1 . This was 
first shown by Jansen and Mazur, who used an exact quantum expansion 
of a in powers of R -1 , 18 This expression will be derived, and 
discussed in some detail in Chapter II. 

The estimate of B £ by Jansen and Mazur makes use of the Unsold 
approximation. 16 In essence this represents an approximation to the 
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resolvent operator, R 0 , in the perturbation expansion, Eq. (l„2l). 
This will also be discussed in more detail in Chapter II. 

In the calculations of Jansen and Mazur, results were obtained 
for both helium and argon, 10 while Buckingham and Pople applied their 
method to more complicated molecules, including those in which a 
dipole moment in one molecule is induced by the permanent quadrupole 
moment of another molecule, 19 > 21 

Since the rare gas experiments had not been performed in 1955, 
the thrust of the theoretical work was to establish experimental 
feasibility. It must be remembered, however, that both of these 
early calculations used approximations to the polarizability which 
were meaningful in the long-range region only, thus the physical 
content of their theories was not complete. 

Additional work on the classical statistical mechanical aspects 
of this problem was done by Hill in 1958, 42 and Ishihara and Hanks in 
1962. 43 

Theoretical interest in the second dielectric virial coefficient 
was renewed by the accurate room temperature experiments of Orcutt 
and Cole in 1967. In particular, the reported sign change (from 
negative to positive) in B e as one goes through the series, (He, 

Ne, Ar) aroused considerable interest, because the early calculations 
predicted positive values for the entire series. This is understandable 
because a negative value for helium means that for some range of 
finite separations the polarizability of the dimer is sufficiently less 



than the sum of the separated atom polarizabilities so that the 
integral which defines is negative. Thus since the early 
calculations used models which at best are valid approximations in 
the dispersion region only* effects in the overlap region could 
outweigh the positive long-range contribution to produce a negative 

value for B . 

e 

Levine and McQuarrie in 1968 modeled the rare gas interaction 
using metallic spheres; 22 DuPre md McTague in 1969 used the 3 £* 
state of H 2 to represent a pair of interacting rare gas atoms. 23 
Both of these calculations indicated that short-range effects did 
produce changes relative to the long-range extrapolations in the 
polarizability curves, but neither result was directly applicable 
to the helium problem in particular. 

Lim, Linder, and Kromhout in 1970 performed both coupled and 
uncoupled Hartree-Fock calculations on the polarizability of the 
He 2 dimer over the range of R from 1.0 to 6.5 a 0 , 24 Their 

failure to produce a negative value for B resulted from the 

e 

failure of their short-range (negative) contribution to compensate 
for the positive contribution in the dispersion region. This gave 
rise to the possibility that previous estimates of the long-range 
effects had been too large. 

For calculations in the low temperature regime, the expression 
for B £ which includes the effects of quantum statistics has been 
recently derived by Bruch, and is given by: 25 
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B e = f dR R2 Aot(R) 9(R) • 0-39) 

Here, g(R) is the quantum-mechanical pair distribution function, 
which reflects the probability that two atoms will be found within 
a separation R of each other. The g(R) is, of course, 
temperature dependent. 

Formal differences between the classical and quantum treatments 
arise because the expression for <M(F)> takes on a somewhat 
different form in quantum mechanics. For this case, <M(F)> is 
expressed as: 

Tr[e* |iH »H(t,F).f)] 

< H (F)> -~j 5 — - ■ (1.40) 

Tr(e sa ) 


where e ■ 1/kT , 


Tr(A) = l <*JA|^> , 

9 * 


and i denotes a possible state of the quantum mechanical system. 
Since H is now an operator, differentiation is accomplished via 
the formula: 


_L e-f 311 _ 

3F e 



j cjH 3H 

d0 e 3? e ] 


(1.41) 
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Here H can be taken to be the Hamiltonian operator for the 
relative motion of a pair of atoms. 

This summarizes the theoretical and experimental work on the 
dielectric virial coefficients up to the point at which the present 
work was begun. 

One goal of the work reported in this thesis is to establish a 
definitive result for the polarizability in the long-range region. 

A second goal is to compute a more accurate value for the short-range 
effect via attainment of the Hartree-Fock limit, in that region, and 
to use these results to compute S £ for both low and high 
temperatures . 

During the time in which this work was in progress, results of 
additional polarizability calculations for He 2 have been reported by 
Buckingham and Watts , 26 and O'Brien et aK 27 Only the latter 
calculation resulted in a negative value for B £ at room temperature. 
No attempt to match the low temperature results of Kerr and Sherman 
was made in either of these calculations. 

D. The Kerr Effect 

In the absence of external fields, the index of refraction 
of gaseous helium is isotropic. However, if one uses polarized light 
to measure the index of refraction, n , when the sample is placed 
in a static electric field, F , he finds that the value of n for 
the case in which the polarization vector of the light and the 
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applied field are parallel (n^) differs from that in which they are 
perpendicular (n^) (field induced birefringence). This effect was 
discovered in glass by Kerr in 1875, and is known as the Kerr 
electro-optic effect. 28 For this phenomenon, the quantity of 
interest is the molecular Kerr constant, m^ which is defined as: 29 


m 


K 


lim 
F + 0 


2(n, - r. x )V m 
27F2 


(1.42) 


In this equation, again is the molar volume. More recently, 
Buckingham and Dunmur have measured m^ for a number of gases. 32 

In order to experimentally determine m^ , one measures the 
field-induced birefringence for a given pressure and temperature, as 
a function of electric field strength. Extrapolation of the ratio, 
2( n H - ^)V m /27F 2 to zero field then yields the value for ny . 

Not surprisingly, it is found that the molecular Kerr constant, 
like the dielectric constant, exhibits a pressure dependence. For 
gases at low density, Buckingham and Pople related m K to molecular 
polarizabilities. 30 Buckingham, 31 and later Buckingham and Dunmur 32 
discussed the corresponding treatment for gases at higher densities, 
by making a virial expansion similar to that used in the dielectric 
constant analysis. In this last paper, it was shown that for light 
of low frequency (much less than the first excitation energy of the 
atoms), the first two terms in the density expansion of m,, for a 

N 

gas of atoms are: 
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m^ ~ p » (1 .43) 

where 

A k = , (1.44) 

and 

B K * K 2 <S<K )> 2 e' U(R>/kT • (!.«) 

Here the p is the density, g(R) is the anisotropy in the 
dipole polarizability for a pair of interacting atoms at a separation 
R , and y is the first nonvanishing hyperpolarizability of an 
isolated atom. This is a high temperature result, because the 
classical limit of the pair distribution function is used. 

For the quantum case, Bruch has shown that replacing the factor 
e U(R)/kT W1 - t j 1 g ( R ) yields an upper bound to the value cf P K , 25 
Differences between the quantum mechanical and classical 
mechanical treatments of the Kerr effect arise because of the 
different form for the ensemble average in the two realms, as was 
observed -in the corresponding treatments of the dielectric virial 
coefficient. The situation is even more complicated in this case 
than it was in the dielectric virial case, because of the need to 
take second derivatives, by applying Eq. (1.41) a second time. 
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Experimental values for and for several rare gases 
at room temperature are given in reference (32), and are listed in 
Table 1.2. Results for for helium gas are conspicuously absent 
from this table. This is because the pressure effect on m„ at 
room temperature is too small to be detected. Since low temperature 
experiments have not been performed as yet, it is not clear whether 
or not a decrease in temperature might render the effect observable. 
One purpose of this thesis is to estimate the magnitude of m v at 
low temperatures. 

There are two additional points regarding the above equations 
worth mentioning: 

1. The field strengths used in these experiments are larger 
than those used in the dielectric constant experiments, 
hence nonlinear polarization effects (via y ) are included. 

2. In the derivation of Eqs. (1 .43)-(l .45) the density 
dependence of y was ignored, so that the higher virial 
coefficients do not involve an integral over y(R) , in 
this approximation. 

Note that the expressions for B (Eq. (1.34)) and B„ are 
similar, except that the respective integrands involve different 
functions of the polarizability components ctg and . Thus, 
regarding the polarizability tensor , the dielectric virial and the 
Kerr virial measurements are complementary. 
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Table 1.2. Experimental first and second Kerr vlrial 
coefficients for several rare gases. 9 


System 

T(°K) 

10 i4 A K (esu) 

' r - --- LU 

10 12 B^(esu) 

He 

296 

0.25 ± 0.02 

— 

Ne 

?96 

0.48 ± 0.04 

— 

Ar 

296 

5.5 ± 0.4 

4.1 ± 0,6 

Kr 

296 

13.0 ± 1.0 

16.0 ± 14.0 

Xe 

296 

36.0 ± 1.0 

65.0 ± 22.0 


a Results are those of reference (2). 
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E. Collision Induced Raman Scattering 

The change in the polarizability during a collision of a pair 
of rare gas atoms gives rise to a light scattering spectrum, which 
is called the collision induced Raman spectrum. 

Experimental detection of such scattering was reported by 
Crawford, Welsh, and Locke 33 as early as 1949 for 0 2 and N 2 gases. 
McTague and Birnbaum 34 in 1968, and Gersten, SI usher, and Surko 35 
in 1970 reported similar observation^ in argon and krypton. Pike 
and Vaughan 36 have recently measured the light scattering spectrum 
of both liquid and gaseous helium, the latter at 4.2°K. 

The expression which relates the Intensity of the scattered 
light, I{u) , to the experimental parameters, has been derived and 
discussed by several workers, and is given by: 3 > 34 > 37 > 81 

I(oj) = F(t)e' itut dt . (1.46) 

—CO 

Here, 

K = V 2 p 2 k 4 1 .7(32*2 e 2 R 2) f (T.47) 

where, V is the scattering volume, p the number density of the 
sample, k^ the propagation vector of the incident beam, 1.. its 
intensity, Ep is the permittivity of the sample, and R 0 is the 
distance from the detector to an arbitrary point in the sample 



(assumed to be much greater than the sample dimensions). The 
expression for F(t) is, for light scattered at right angles. 
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F(t) 


= <Aa(t=0)*Aa(t)>(n i *n f ) 2 

0 . 4 $) 

+ <e(t=O)6(t)( 1 / 4(3cos 2 0(t)-l))>{45 _1 (3 + (n i ‘n f ) 2 )} . 


In this equation, n. and are unit polarization vectors of the 
incident and scattered beams respectively; 0 ( t ) is the angle 
between the internuclear vector of a diatom at time t and an axis 

defined by its vector at time t-0 ; Aa(t) and 3 ( t ) are 

respectively the mean incremental polarizability and anisotropy of 
the diatom at time t , while the brackets denote an ensemble average. 

Experimentally, one often measures depolarization ratios which 
may be defined in a variety of ways, which depend on the geometry of 
the experiment. In the low temperature work of Pike and Vaughan, 
light from an argon ion laser was passed horizontally through a cell 
containing gaseous helium. 36 The incident light was introduced in 
either horizontal (H) or vertical (V) polarizations, and the 
Intensity of the light scattered at 90° was analysed in both the H 
and V polarizations. Their data for helium gas at 0.96 atm and 

4.?2K is given in Table 1.3 for a frequency shift of 21 cm’ 1 from 

thq exciting line, using a spectrometer with an instrumental 
linewidth of 5 cm" 1 . 
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Table 1.3. Experimental depolarization ratios at a shift of 30K in 
helium liquid and gas. 3,15 


Depolarization Ratio 

Liquid 

Gas 

V-V/V-H 

1.52 ± 0.2 

1.33 ± 0.2 

H-V/H-H 

1.1 ±0.2 

1.0 ± 0.2 

V-V/H-V 

1.29 ± 0.1 

1.39 ± 0.2 

a Reference (36). 



b Gas experiments were 

done at 4.22°K at a pressure 

of 0.96 atm. 


Liquid experiments were done at 4.22°K and 1.15 atm. 
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Eqs. (1.46) and (1.48) give the following expressions for the 
depolarization ratios measured by Pike and Vaughan. 



V-V/V-H 

- [A(w) + 


0] . 

(1.49) 


H-V/H-H 

= 1 , 



(1.50) 


V-V/H-V 

= V-V/V-H 

4 


(1.51) 

Here, 

the Fourier transforms A( w ) 

and B(u) are 

defined by: 



A(o>) 

= |°° <Aa{0)Aot(t)>e' iwt dt 

— oo 

» 

(1.52) 

and 







B(») = | 

00 

<e(o)s(t)p 

2 (cose(t))>e" ltot 

dt . 

(1.53) 


Here, P 2 U) - /£(3x 2 -l) is the second Legendre polynomial. 

We see that the data in Table 1.3 are consistent with these formulas, 
and, furthermore, imply that A(a>) << B (to) since in this limit 


V-V/V-H = % . 


(1.54) 
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Although the accuracy of the reported data is not sufficient to 
distinguish the contribution of A(w) to the depolarization ratios, 
we shall estimate its effect computationally in order to verify that 
it is not observcble in the present experiments and to establish its 
magnitude for use in planning future experiments. To avoid dealing 
with time-dependent correlation functions, which are difficult to 
compute, we shall instead deal with integrated intensities, or the 
zeroth frequency moments of Eqs. (1.52) and (1,53) 


A 



A(w) dto = <Aa(0) 2 > 


(1.55) 


and 


Eq. 


B 



B(w) 


du 


<B(0) 2 > . 


(1.56) 


'o 

We shall calculate the depolarization ratio corresponding to 
(1.49) and defined by 


[A + VsBl/CVsB] , 

(1.57) 

% + 15A/B . 

(1.58) 


Thus the interesting part of D involves averages over both the 
incremental mean polarizability Aa and the anisotropy 6 . Higher 
moments of (1.52) and (1.53) lead to other experimental parameters 



34 


such as linewidths, and approximations to line shapes. These 
extensions, however, involve more than the simple integrals of Eqs. 
(1.55) and (1.56). 


Semi -empirical work by Levine and Birrbaum on the collision 
induced Raman scattering of rare gases has been directed toward 
obtaining a model for e(R) . Using a modification of the point 
dipole anisotropy, 37 

e(R) = 6ag R" 3 + BR' P , (1.59) 

they derived high temperature, two-body expressions for the first two 
moments of the spectral functions, A(u) and B(o>) , and obtained 
values for p and B by comparison to the experimental results for 
these moments, for argon, krypton, and xenon. Typical values so 
obtained for argon are on the order of 10 and -1.0 x lo 5 A 13 
respectively. This is then interpreted to mean that short-range 
effects as represented by the Br -p term are indeed significant. 

Recently, Gelbart has given a quantum mechanical treatment of the 
depolarization and inelastic scattering of light by gases and liquids 
composed of atoms and/or isotropically polarizable molecules. 38 
Effects due to collisions of more than two particles were also included. 

F. Perspective 

Apart from stricly formal treatises, prior theoretical work on 
estimating the properties of interest has been restricted to the 
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high temperature regime, and have involved the coupling of 
polarizabilities computed at different levels of approximation with 
classical averaging procedures. The various polarizability models 
which were used and the conclusions of their respective calculations 
are discussed in some detail in subsequent sections; however, for 
purposes of understanding the motivation for the present work, we 
note at this point that the most successful of these model approaches 
employed the Hartree-Fock approximation. 

The meaningful use of any model depends on the existence of 
well-defined, reproducible model limits. In the case of the 
Hartree-Fock approximation, this requires assurances that in a given 
calculation the Hartree-Fock limit (the limit of zero basis set 
error in the Hartree-Fock-Roothaan scheme) has been reached . 39 This 
is a particularly important point when one considers the history of 
He 2 pair potential calculations. 

In an early Hartree-Fock calculation on the ground state 
potential energy curve of He 2 , Ransil reported that the dispersion 
minimum was observed . 40 This minimum was later shown, by further 
calculations, to be spurious, and was attributed to an idiosyncracy 
in the basis set which Ransil used . 41 Thus the validity of any 
conclusions drawn from the use of the Hartree-Fock-Roothaan scheme 
depends directly on the closeness of the calculation to the 
Hartree-Fock limit. 
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None of the Hartree-Fock calculations of the polarizability 
which have been thus far reported, purport to represent the 
Hartree-Fock limit . 24 > 26 > 27 It was our intent in carrying out the 
present work to remove this uncertainty and establish the need or 
lack thereof of more accurate calculations, which include electron 
correlation effects. 

It goes without saying that even if a certain model duplicates 
experiment at one temperature, it does not necessarily follow that 
the same model will be equally useful at another temperature for 
which the physics of the problem might be different. This means 
that a given polarizability result must be tested at both high 
temperatures using classical statistical mechanics and low temperatures 
using quantum statistical mechanics if one wishes to make claims 
regarding the usefulness of the model from which the polarizability 
resulted. In addition, isotopic effects may be introduced in the low 
temperature region so that the ability of a given model to produce 
such effects must also be guaranteed. 

Since no low temperature work had been carried out up to now, 
previous workers did not subject their polarizability results to 
these important tests. This is, of course, insignificant when 
considering models which make no claim for general usefulness. 

The quantum statistical methods seem to be on somewhat firmer 
ground than the polarizability calculations, although they have 
never been used for anything other than equation of state 
applications. 
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It was the intent of the work reported in this thesis to couple 
the Hartree-Fock limit polarizabilities to quantum statistical 
mechanical results for the two helium isotopes ( 3 He and 4 He) in 
order to duplicate the results of both low and high temperature 
experiments. This was done with the intent of establishing a basis 
on which one can judge whether or not he is justified in using the 
same approximations in the calculation of the polarizabilities of 
other systems or in the calculation of different but similar 
molecular properties. 



CHAPTER II 


POLARIZABILITIES OF DIATOMIC MOLECULES IN THE LIMIT 
OF LARGE INTERNUCLEAR SEPARATIONS 

Early estimates of the dielectric virial coefficient employed 
extrapolations of approximate polarizabilities from the dispersion 
region to the regions of smaller R . 18,19 Q ne such approximation, 
the point dipole model, has already been discussed. It is apparent 
that if the leading contribution to the point dipole incremental 
mean polarizability (see Eqs. (1.37) and (1.38)), 

MR) = 4c*q R -6 + . . . , (2J) 

is used for all values of R to compute a value for B £ for a gas 
of ground state atoms, a positive result is obtained, in contrast to 
the observed negative result for helium. It is natural to ask 
whether this discrepancy is due to the inaccuracy of the point dipole 
approximation itself, or due to the use of a long-range result for 
all values of R . 

In this chapter, we show that the point dipole model gives the 
leading asymptotic behavior of Aa , although the magnitude of the 
leading multipole coefficient is too small by about a factor of four. 

38 
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We derive the long-range multipole expansion of the diatomic 
polarizability in Section A, illustrate its connection to the point 
dipole model in Section B, and discuss the calculation of the 
leading terms for the H 2 , HeH, and He 2 systems in Section C. The 
results of the first two sections were first obtained by Jansen 
and Mazur. 18 

A. The Long-Range Multipole Expansion of a(R) 

A rigorous treatment of the pair polarizability of two atoms in 
the long-range region can be carried out using Rayleigh-Schrodinger 
perturbation theory. In the long-range region, the electron clouds 
of the two atoms do not overlap, and it is not necessary to require 
the diatomic wavefunction to satisfy the Pauli Principle for electrons 
centered on different atoms. This greatly simplifies the calculation 
of the polarizability. 

Let us consider two S-state atoms, a and b , separated by a 
distance R . The starting point of the development is the general 
expression for the polarizability, Eq. (1.21), 

a = -2<¥ 0 |y R 0 u|V . (2.2) 

In this expression, Vq is the wavefunction for the diatom in the 
absence of the electric field, Rq is the reduced resolvent for the 
diatom, symbolically. 
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Ro = (E 0 - Ho)' 1 , (2.3) 

where H 0 is the diatomic Hamiltonian and E 0 is the ground state 
energy s 

(H 0 - E 0 )f 0 - 0 . (2.4) 

The u is the total dipole moment operator, which is a sum of the 
dipole operators for atoms a and b , 

y " y a + y* 3 . (2.5) 

In the following derivation, we shall consider separately the cases 
that the external field is along the internuclear axis (designated 
the z-axis) or perpendicular to it. For simplicity of notation, the 
expansions of Eqs. (1.18) and (1.21) shall be written as scalar 
expansions with the understanding that the resulting polarizability 
expressions apply for the respective components. 

We now wish to obtain an expansion for a in inverse powers of 
the separation R . The R-dependence in a occurs both in the 
wavefunction and the resolvent Rq . Vie proceed by subdividing 
the diatomic Hamiltonian according to 

Ho 


“ h 0 + V , 


( 2 . 6 ) 



where ho is the sum of the atomic Hamiltonians for a and b , 
and V represents the Coulombic interaction between a and b . 
Similarly, for the energy, we have 


E 0 - e 0 + AE . 


(2 = 7) 

To obtain the multipole expansion of Tq » we first resolve the 
Schrodinger Equation (2,4) into perturbation equations based on 
Eq. (2.6), 

*0 = <f>0 + 4 ]) + 4^ + ... 

* 5 

(2.8) 

where 



AE = E^ ^ + E^ + .o. , 


(2.9) 

( h Q - C Q ) 0 ~ 0 , 


(2.10) 

(ho - e 0 )4 1} + (V ‘ E (1) Ho - 

0 , 

(2.11) 

(h 0 - e 0 }4 2 ) + (V - E^ 1 

<f>o “ 0 > 

(2.12) 


with the full normalization conditions, 
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<$oK (1)> = 0 , (2.13) 

<+ol*o (2> > “ <?o ( 2 ) l*o> « -k < 4 1 ) l* 0 1)> • (2.14) 

The solution of these perturbation equations is well-known. 

The <p 0 is a simple product of atomic wavefunctions for atoms a 
and b , 

♦ o = a o t>o • (2.15) 

Introducing the reduced resolvent for h 0 , 

ro = ( e o ‘ ^ o > ~ 1 3 (2.16) 

we have 

55 r o(V - E (1 ^ )(j) 0 , (2.17) 

and 

¥ o < 2> * r ° (V " e( 1 ) ) ¥ o ( 1) ~ H*o<^ 1 ) l^ 1)> . (2.18) 

The R-dependence is still implicit in these equations. However, 
both V and AE have well-known multipole expansions in inverse 
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powers of R . 82 For our purposes we need consider only the leading 
terms , 


V o 

\f - + 

* 1 * ■ * 9 


and 


= 0 , 


•( 2 ) . 


. bL + 

R6 * • * » 


where V 3 is the dipole-dipole interaction operator. 


1 / _ a 3 A b , A a -b -b 

Va - “x “x + •V v y - v z • 


(2.19) 


(2.20) 

( 2 . 21 ) 


( 2 . 22 ) 


and 


Ce = -<i>o|V 3 r 0 V 3 |<f>o> . (2.23) 

Then, introducing the multipole expansions into Eqs. (2.17) and 
(2.18) and collecting terms, 

= 4>o + R" 3 r 0 V 3 (j) 0 + R' 6 [r 0 V 3 r 0 V 3 - ^<+olV 3 roV 3 |* 0 >]*o + 


(2.24) 
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Although we have not exhibited all the terms in f 0 , we have shown 
those terms which contribute to a through order R“ 6 . 

Finally, the R-dependence of the resolvent R 0 is obtained by 
invoking the operator identity, 

(A-B) -1 = A " 1 + A" 1 B (A-B ) ~ 1 . (2,25) 

for A = eo - ho and B = V - aE . Thus, 

Ro = r Q + r 0 (V - AE)R 0 . (2.26) 

Iterating this equation we obtain 


Ro = r 0 + r 0 (V - AE)r 0 4 r 0 (V - AE)r 0 (V - AE)r 0 + , 

(2.27) 

and then introducing the multi pole expansion, 

R o = r 0 + R" 3 r 0 V 3 r 0 + R“ 6 [-C 6 rQ + r 0 V 3 r 0 V 3 r 0 ] + .... (2.28) 


Again we have exhibited only those terms which contribute to a 
through order R' 6 . 

The multipole expansion of a may now be obtained by 
substituting the expansion for y 0 » Eq. (2.24.), and for R 0 , Eq. 
(2.28) into Eq. (2.2). When this is done, we obtain 
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a 


= 2a 0 + pg + 



(2.29) 


where 


2a o = -2<0 q|p r 0 y|i> 0 > » (2.30) 

A 3 - - 4< <i> o I V 3 r 0 y r 0 y|<j>o > - 2<<f>o|y r 0 V 3 r 0 y 1 4>o> » 

(2.31) 


end 

A 6 = -4<<j>o |V 3 r 0 V 3 r 0 y rg y^o 1 * “ ^ < 4>o [ ^ 3 r 0 y r 0 V 3 r 0 y 1 4>o > 


- 2<<|>o jy r 0 V 3 r 0 V 3 r 0 y 1 4>o> ' 2«t> 0 |V 3 r 0 v r 0 y r 0 V 3 j<j> 0 > 

2 A 
+ 2 <vf) 0 j V ^ r 0 V 3 I 0 ><t t‘O 1 r 0 Vi|<j>0 > 

+ 2 <<J>o|V 3 r 0 V 3 j 4> 0 > < <> 0 |y r o ^ 1 4 ) o > » 

(2.32) 


where these expressions apply for the parallel component of a 
when y - y v and the perpendicular when y = y . 

*— A 

We note that the individual components of a involve terms 
which decrease as R -3 . In the next section, however, we show that 
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A, 



(2.33) 


Thus, i-he incremental mean polarizability has the asymptotic behavior. 


Aa 



2a 0 


+ 


aA 6 
R6 + 


(2.34) 


where 


aA 6 = 



(2.35) 


The terms of a which decrease as R" 3 yield the leading terms in 
the asymptotic expansion of the anisotropy 


0 


01 II “ a l 






(2.36) 


For reference in later sections of this chapter, we note that 
the reduced resolvent operator, r^ , can be expanded in terms of 
the excited eigenstates of h 0 = hj + hp : 


r 0 = J' l a K b 0 ><a K b 0 1 ] y, 

K e a -e a t 

K 0 L 


a 0 b L ><a 0 b L | + y , y , J a K b L ><a K b Ll 
H L L a a . k "k 


£ Nn 

L 0 


K L 


(2.37) 


where a K and b L represent excited states on the respective centers 
and the prime means that the terms K = L = 0 are omitted from the sums 
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B . Co nnection between the Point Dipole Model and the Quantum 
Mechanical Multipole Expansion of the Polarizability 
It is instructive to relate the exact long-range behavior of 
the polarizability to the predictions of the point dipole model. 

This is the object of the present section. 

Before treating the problem mathematically, we can make certain 
qualitative statements regarding what to expect. Since the point 
dipole model ignores the details of atomic structure, whereas the 
quantum mechanical treatment correlates interatomic electronic motions, 
we shoulo expect differences in those terms which are present in both 
the field-free and field-dependent quantum multipole expansions of 
the interaction energy. On the other hand, since the point dipole 
model reflects modifications in the local field at a given atom due 
to the field-induced dipole moment in the other atom we might expect 
some degree of similarity between the terms of the two expansions 
which do not appear in the field-free quantum expansion, and represent 
therefore field-induced terms in the quantum expression. 

The manifestation of the first of these points will be clearly 
seen when we consider the actual results of the calculation of the 
long-range coefficients for specific systems. The second point is 
immediately accessible, and will now be illustrated by showing that 
the A 3 term which results from the point dipole model is identical 
to that which results from the exact quantum mechanical multipole 
expansion. For simplicity, we consider the aJ term only. Since 
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the components of the dipole moment operator behave like the 
corresponding components of the vector r , and since y 2 = y z + y^ , 


we can write Eq. (2.31) as : 

A 3 = -2[-4<^o|u* y z r 0 y 2 r 0 y 2 j<f> 0 > - 2<0 o |y z r 0 y^ pjj r ° 


= [Tl v T2] . 


(2.38) 


Further, only those parts of r 0 which are of p-type symmetry (behave 

J 

like Y 10 spherical harmonics) can lead to nonzero contributions to 
A 3 . Thus , using superscripts to label the respective centers , 
subscripts to label different atomic states of a given symmetry, 

Ae^ to denote e 0 - , and writing <)> 0 = a 0 b 0 , we have: 


Tl 


sp 

K,L,M 


+ 8 p 

K,L,M 


<a o^o 1 PzUz 1 PrPl >< Pk1 1 a o ><: PL I Pm >< PmI vt\b 0 > 

(Ae| + Aej?) &e b 


(2.39) 


<a obp l^zn 2 |PKPL ><: PLln 2 jb 0 ><p|| p§><p§|y||a c 
(Ae| + Ae^Ae* 


Because of the orthonormality of the atomic states, this last 
equation can be factored, using the identity 


(ab) -1 = [{a+b)a] _1 + [(a+b)b] -1 , 
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T1 - o v i r <a Q I iPK > 2<b O 1 Vz I Pl > 2 + <a 0 1 Wz | PK >2<b o | £z | PL> 2 -, 


1 M v <a ? |il ! pg>2 ][P *■!»] , 

K Ae K L tei. 


~b i _b 


(2.41) 


O a b 

~ 2ct o a 0 • 


(2.42) 


Similarly, 


T2 


<a Q |£flPK > 2 < PLUzl b o >2 


4 l Q I H z l TK' ^Li^zl^O' + 4 j 


K,L 


4e k 


<b o I U^l P l» z ^PkI S z l a 0 >2 


K,L 


Ae L Ae K 


(2.42) 


o ^ b 

2 «o °o • 


(2.43) 


Thus , 


A3 = 4ag oq , (2.44) 

which is identical to the general point dipole result for a 
hpterppolar diatom. 

In a similar fashion, using the symmetry reduction procedure 
discussed in tyie next section and described in Appendix A, one can 
show that. 
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a£ = -2a a 0 aj . (2.45) 

Thus, as was first shown by Jansen and Mazur, the leading term in 
the point dipole expansion of the polarizability is identical to 
that obtained from the exact quantum mechanical multipole expansion. 10 
A significant implication of this is that the leading contribution to 
the mean incremental polarizability for a diatom goes as R" 6 , since 
A3 + 2A^ = 0 . 

It does not appear possible to obtain a simple relation between 
the quantum mechanical and point dipole expressions for the higher 
multi pole coefficients, however. The accurate numerical evaluation 
of the terms A$ and Ag is described in the next section. 

C. Accurate Long-Range Polarizability Calculations for the 
H 2 , HeH, and He?. Systems 

The initial application of the equations of Section A to the 
calculation of the long-range contribution to diatomic polarizabilities 
were made by Jansen and Mazur in 1955. 18 In particular, they used the 
Unsold approximation to simplify the use of Eq. (2.32), and computed a 
value of the second dielectric virial coefficient for helium gas. The 
Unsold approximation involves the replacement of the energy 
denominators in ro (see Eq. (2.37)) with some average value, Ae , 
so that r 0 can be written: 



51 


ro 



(2.46) 


and using the completeness of the states ip K , 

r 0 = — (1 - | <f>o ><( f>o I ) • (2.47) 

Ae 

With this approximation, one avoids having to compute the spectrum 
of atomic excited states. 

The results of the Unsold approximation calculations for helium 
yield values of Ag and Ag of 7.59 and 3.94 oq respectively; 
here, a 0 is the helium atom polarizability. 

Thus, as far as approximate calculations in the dispersion region 
are concerned both the point dipole and the Unsold approximation yield 
positive values for for helium gas. However the possibility of 
obtaining a negative value via the extrapolation of an accurate 
dispersion polarizability has not been completely eliminated. To 
test this, we have computed A 6 as given by Eq. (2.32). 

The only accurate calculation of A 6 prior to the present one 
was that of Tulub ejt al_. for H 2 , in 1970. 44 In their work, they used 
a multipole expansion of the interaction energy of a pair of hydrogen 
atoms in an electric field, and averaged this energy over all 
orientations of the internuclear axis with respect to the field. 

They then solved the resulting perturbation equations variationally. 

A drawback of this procedure is that only the mean value of A 6 is 
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obtained; hence strictly speaking, no accurate detailed knowledge of 
the A 6 tensor was available at the time the present work was 
performed, even for a system as simple as H 2 ! 

As in the Unsold approximation, Eq, {2.32) serves as the basis 
for an accurate calculation in the long-range region. Ignoring, for 
the moment the question of how r 0 is to be obtained, we note that 
Eq. (2.32) is not in a computationally convenient form because the 
various terms contain all the Cartesian components of y in V 3 . 
Eq. (2.32), however, can be simplified so as to yield an expression 
which involves only the z-components of the various operators. This 
symmetry reduction consists of a straightforward but tedious 
application of the Wigner-Eckardt theorem^ 5 to the terms of the 
equation which results when Eqs. (2.22) and (2.37) are introduced 
into Eq. (2.32). A discussion of this simplification is given in 
Appendix A. 

The expression for the components of A 6 which results can be 
wri tten : 

Ac = -(1 + P ab )[Mi<Vr 0 Vr 0 y r 0 y> g + M 2 <Vr 0 Vr 0 y r 0 y> D 

+ M 3 <Vr 0 y r 0 V r 0 y> s + My<Vr 0 y r 0 V r 0 y> D 
+ M 5 <y r 0 Vr 0 Vr 0 y> g + M 6 <y r 0 Vr 0 Vr 0 y> D (2.48) 

+ M 7 <Vr 0 y r 0 y r 0 V> s + M a <Vr 0 y r 0 y r 0 V> D 
+ M 9 <y r 0 y><Vr§ V> + M 10 <y r a y><Vr 0 V>] . 
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In this equation, V - u 3 , u - u 3 , P is the permutation 
operator for the two centers, and the brackets denote an expectation 
value with respect to <j> 0 , The coefficients are given in 
Table II. 1. The subscripts S and D are explained below. 

There are several points worth noting about Eq. (2.48). First, 
its use requires only atomic information as input. That is, r 0 
refers to the reduced resolvent operator for the Hamiltonian, 
ho = h 3 + bo » s0 that the excited states which enter the spectral 
expansion of r 0 , Eq. (2.37), a^e products of atomic eigenstates on 
the respective centers. This is a significant simplification as far 
as computations are concerned, because atomic calculations are much 
easier to carry out than the corresponding calculations on the 
diatomic systems. Second, because of the orthonormality of the 
atomic states, only certain symmetries enter a given r 0 so as to 
yield nonzero contributions to A 6 < The particular symmetries 
required in a given r 0 depend on its position in the matrix 
element under consideration. As an example, consider the first term 
in the above expansion, and assume that the resolvents are numbered 
from left to right. Since V - p 3 p^ , the only nonzero contributions 
from the first Tq are due to terms which denote states of p-symmetry 
on each center. Similarly, the second resolvent can only involve 
those states which are represented by products of excited S or D 
functions on each center. Lastly, the third resolvent involves only 
those products which consist of p-states on center a and b 0 on 
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center b . This last statement is appropriately modified under the 
action of P ab by using excited p-states on center b , and 3 q 
on center a . 

The above points allow for a notational simplification which 
unambiguously defines the states entering each resolvent. This 
consists of explicitly specifying the symmetry (S or D) of the middle 
resolvent, so that the symmetry of the other two can be obtained via 
the orthogonality of the atomic states. Thus, the first two terms 
of Eq. (2.48) result from the two possibilities for the second r 0 
which enter a single matrix element, where the subscript refers to 
the symmetry of the middle resolvent. 

The computational ease with which Eq. (2.48) can be used rests 
on one's ability to compute the atomic resolvents. In most cases, 
the exact calculation of these operators is, of course, impossible. 
Thus, one generally resorts to approximate schemes to compute them. 
The method used in this work was to obtain "variational" 
approximations to the excited spectra by diagonalizing the matrix of 
the Hamiltonian operator in a basis set, for the atomic systems of 
interest. 

In our work, the individual components Ag and Ag were 
obtained by performing separate calculations for the S , P , and 
D parts of the resolvents. Moreover, since we also symmetry reduced 
the expression for Ag" , only the = 0 components of the atomic 
states were needed. 
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Table II. 1, The coefficients s M n » appearing in the expression for 
Ap and A & L . a 



Mi 

m 2 

m 3 

m 4 

Ms 

m 6 

My 

Me 

B 

Mio 

Al 

24 

12 

16 

22 

8 

11 

8 

n 

-12 

-12 

Ac 

24 

-6 

4 

19 

2 


2 

% 

-12 

-12 


a See Eq, (2.48). 
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The atomic basis functions for the helium calculations were 
taken to be of the form, 


^nJl ~ a ° ’ (2A9) 

where ao is Hart and Herzberg‘s twenty term, directly correlated, 
variational approximation to the He atom ground state wavefunction,, 46 
The form of this wavefunction and the parameters used in it are given 
in Table 11,2. 

The advantage of choosing this particular form for the basis 
set is that, for sufficiently accurate a 0 s the matrix elements of 
the operator (hg - e 0 ) are given, to sufficient accuracy by: 47 ’ 48 

2 

" e o)l* nA > = I <a ol v i F nH , VnJ a O > • ( 2 - 50 ) 

i_l 

This last equation is exact only if a 0 is an exact eigenfunction 
of h^ . Thus, this method is equivalent to replacing h* by hf , 
so that for the particular a 0 used, one has (h® - e 0 ) a o * 0 , 
while ignoring (h® - fij) * The justification for the use of this 
approximation in conjunction with the Hart-Herzberg wavefunction 
will be given later by comparing the results of expectation values 
of various operators for the ground state of He which were computed 
using it, to the corresponding results obtained by using more 


accurate wavefunctions. 
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Table II. 2. Numerical parameters for the Hart-Herzberg He atom 
wavef unction. 3 


m 


m 


m 


m 


m 


m 


1 0,33729424 

2 0.08088340 

3 -0.21312975 

4 0,02003854 

5 -0,02871601 

6 -0,01543812 

7 -9.2189670 x 10“ 3 


8 4.3290467 * 10" 3 

9 7.8697645 x 10" 4 

10 -1.7755868xl0- 3 

11 -7.4084122 x 10“ 4 

12 1.6305837 x 10" 6 

13 -2.7310618 x 10" 4 


14 

6.2740944 x 

TO" 6 

15 

-6.3829315 x 

10" 5 

16 

-1,8442313 x 

10“ 4 

17 

1.5585701 x 

10“5 

18 

6.4835063 x 

10' 4 

19 

6.8860244 x 

10' 4 


a o “ Ne ^ S [l + C]U + c 2 t 2 + C3S + C4S 2 + C5U 2 + CgSU + Cyt 2 U 

+ CgU 3 + c 9 t 2 u 2 + c 10 st 2 + c n s 3 + c 12 t 2 u 4 + c 13 u 4 + c ll4 u 5 
+ C ls t 2 u 3 + c 16 s 2 t 2 + c 17 s 4 + C le t 2 U + c 19 t 4 ] , 
where, N = 0,07164092 , k - 1,35 . The coordinates, s, t, and u 
are defined in terms of r 1$ r 2 , and r 12 by: 


s = r 1 +r 2 , t « r l -r 2 , and u = r 12 . 
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For these calculations, the form of the F nl 's for the helium 
atom was taken from a calculation of the hyperpolarizabili'ty of 
helium, 47 It consists of the following functions: 

S-symmetry: (r^ -2 + r£~ 2 )ja 0 > 2 < n £ 4 

(r?~ 2 + r 2 ~ 2 ) r i r 2 cose 12 |a 0 > 2 < m < 5 

P- symmetry : (r^ -1 z l t r*” 1 z 2 )|a Q > 1 < n £ 6 


D-symmetry: (r?- 2 (z 2 - V 3 r 2 ) + r»- 2 (z 2 - 1 / 3 r|))ja 0 > 2 < n < 4 

( r m -2 + r™ -2 )( z 1 z z - l ^r r r 2 cose 12 )|a 0 > 2 < m <_ 5 

In order to obtain a final basis for the helium atom 
calculations, the convergence of Ag 1 and was monitored, along 
with that for the atomic hyperpolarizability, polarizability, and 
C 6 coefficient. The first of these atomic properties allowed an 
assessment of the completeness of that part of the resolvent with 
S- and D-symmetry. The second and third of these additional 
properties serve to test the completeness of the P part of r 0 . 
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In the latter case, the utility of these specific properties as test 
properties is especially pronounced, because one obtains a rigorous 
upper bound to both the polarizability and the C 6 coefficient 
corresponding to the Hart-Herzberg a 0 , as they each correspond to 
a second order energy. The final basis set was obtained by varying 
the number of functions of each symmetry type until convergence for 
all these quantities was achieved. The convergence of the atomic 
hyperpolarizability, yo > a£ , and Ag as a function of the 
number of terms in the basis set is given in Table II. 3. 

Besides testing the accuracy of our calculation in the above 
manner, sum rule checks were employed to test the internal consistency 
of our calculations. It has been shown by Dalgarno and Epstein 49 that 
if <j) 0 is an exact eigenfunction of h 0 , and if the basis set for 
a variational calculation of the states of h 0 contains the function 
<f>o » then the following sum rule holds: 

1 l<+ol0|*o > l 2 = < *ol0 2 l+o> • (2.51) 

n 

In our case, the Hamiltonian matrix for the linear variational method 
was constructed as if the Hart-Herzberg wavefunction was exact. That 
is, the Hamiltonian matrix is that which refers to h 0 . Hence, if 
our calculation is numerically precise, and internally consistent, 
the above sum rule should hold for operators 0 , which are of the 
form of our F 's . Since the F .'s were different for each 

n£ nJi 
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symmetry, the numerical consistency of each part of the calculation 
could be independently tested. The results for the sum rule checks 
which we performed are given in Table I I. 4. 

In Table II. 5, we present the results of calculations for 
expectation values of various radial operators. The good agreement 
between the results obtained by using the Hart-Herzberg ground state 
wavefunction and those obtained through the use of the very accurate 
(1,078 term) Pekeris wavefunction 50 testifies to the good quality of 
the former, and justifies the method which we used to compute the 
matrix elements of the Hamiltonian. 

All the integrals needed in this calculation were obtained using 
the method of Calais and Lowdin. 51 For the most part, their method 
can be used without modification. However, for certain of the 
integrals needed in the construction of the Hamiltonian matrix, 
modifications had to be made. The analytical evaluation of these 
integrals is discussed in Appendix B. 

For the hydrogen calculations, the basis functions <b „ for 

nt 

the excited states were taken to be of the form, 

+nA " F nll *0 . 



r a Vr> 


» 


where 
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Table 

II. 3. 

Convergence of y 0 , A J , and Ag for He, as a 
function of the number of terms in the basis. 


S Type 

P Type 

D Type 


J 

K 

As 

Ni 

n 2 

N3 

n 4 

n 5 

5 

0 

5 

5 

0 

38.9 

54.98 

25.61 

3 

3 

6 

3 

3 

38.7 

57.71 

26.56 

4 

3 

6 

4 

3 

42.5 

60.09 

28.42 

5 

3 

6 

5 

3 

42,6 

60.09 

28.43 

6 

3 

6 

6 

3 

42.6 

60.09 

28.43 

4 

4 

6 

4 

4 

42.5 

60.09 

28.43 

5 

3 

7 

5 

3 

42.6 

60.09 

28.43 


Table II. 4. Sum-rule check for the operator 0 (see Eq. (2.52)), 


0 

ri+r 2 

2 2 
r l +r 2 

P 2 

z 2 M-H(r z M) 

ZiZ 2 -H Y 'i r 2 cose 12 

ZI<0|o’|n>l 2 

0.32155684 

0.82473581 

0.75239640 

0.62955186 

0.35870841 

<0|0 2 |0> 

0.32155688 

0.82473588 

0.75239638 

0.62955184 

0.35870840 


CD 

ro 



Table II. 5. Calculation of properties of a ground-state helium atom; all properties in atomic units. 

N g is the number of electrons, C 6 is the coefficient of the R" 6 dispersion energy, 
ao is the static polarizability, y 0 is the static hyperpolarizability . 


Property 

C 6 

«0 

TO 

2 2 
< r i +r 2 > 

<r l *r 2 > 

N e 

m 

O 

o 

w 

Value, 
this calc 

1.458 

o% 

r^ 

m 

0 

42.6 

2.3867 

-0.064736 

2.0000000 

-2.9037179 

Accurate 

value 

1.4605 ± 0.0025 a 

1.380, b 
1 . 383 c 

42. 8, d 
42.6 e 
43. l e 

2.3870 f 

-0.064737 f 

2 

-2.9037244 f 


a M. H. Alexander, 0. Chem. Phys. 52, 3354 (1970). 
b W. D. Davidson, Proc. Phys. Soc., (London) 87, 133 (1966). 


c Reference (7). 

d M. W. Grasso, K. T. Chung, and R. P. Hurst, Phys. Rev. T_67_, 1 (1968). 
e Reference (47). 

f 

Reference (50). 



Here t|> 0 is, of course, the exact ground state wavefunction. It 
was found that n £6 for each symmetry was sufficient to obtain 
convergence for Ag and Ag . 

For the helium- hydrogen problem, the final basis sets for the 
above respective parts were used together. 

The values of A l and Ag obtained in this calculation for 
the H 2 , HcH, and He 2 systems are given in Table H.6. 

After this work was completed, Buckingham et al_. reported the 
results of a variational calculation of both A 6 and A 0 for H 2 . 52 
This calculation was virtually identical to ours in method, and 
their A 0 results are in excellent agreement with ours. For 
reference, we include these results, those of reference (44), the 
point dipole, and the Unsold results in Table II. 6. 

If the results of Table I1.6are used to compute a value of B £ 

for helium, a positive value again results. Thus, one concludes 

that calculations which take only long-range interactions into 

account are not sufficient. This means that, not surprisingly, the 

overlap region contributes significantly to the measured value of 

B £ , so that ot(R) must be computed over the entire range of R if 

a meaningful value of B is to result. 

£ 
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Table 1 1. 6. Accurate 

and approxmate values 

£ A 11 

of A g and 

A^ for H 2 , 

HeH, and 

He 2 j units of A 6 are 

9 , . 

a 0 , units 

of Ag are a 0 . 


h 2 

HeH 

He 2 

Ag(accurate) 3 

2558. 59 e 
(2558.59) c ’ e 

517,472 

60.0922 

Ag(accurate) a 

1268, 25 e 
(1268.25)c s e 

258.871 

28.4308 

Ag (point dipole) 

729.000 

146.451 

21,1620 

Appoint dipole) 

182.250 

36,6129 

5.29050 

Ag(Unsold) b 

1030.85 


20.0874 

Ag(Unsold) b 

552.445 


10.4157 

a' Ic 

90639.5 



A, 1 ' 

22010.3 




The number of figures reported in this table corresponding to the 
accurate calculation (see Eq. (2.48)) reflect the numerical precision 
obtained. Except for the H 2 calculations, the error introduced 
through the use of h 0 (see Eq. (2.50)) results in about four 
significant figures of accuracy. 
b Reference (18). 
c Reference (52). 

e Reference (44) reports A = (A 6 + 2A 6 )/3 to be 1698 a^ . 


CHAPTER III 


POLARIZABILITIES OF DIATOMIC MOLECULES 
IN THE OVERLAP REGION 


From the results of the previous chapter, It is clear that 
consideration of overlap effects is necessary if one is to accurately 
calculate the incremental polarizability of two helium atoms. We 
shall somewhat arbitrarily define the overlap region for helium to 
span separations from 4a 0 to 8a 0 , which covers most of the van der 
Waals well region. In the overlap region, we shall apply the 
Hartree-Fock model to calculate the polarizabil ity tensor of He 2 at 
four internuclear separations: 4, 4.7, 5.5, and 6.0 a 0 . We shall 
also show that the Hartree-Fock model reduces to the point dipole 
model at large internuclear separations, where orbital overlap is no 
longer Important. Our results are presented in Sections C-F. 

Before discussing our own work, we review previous calculations 
of the helium pair polarizability in the next section, and review 
coupled Hartree-Fock theory in Section B. 

A. Previous Calculations in the Overlap Region 

Early attempts at estimating the effects of overlap on the 
polarizabilities of rare-gas diatoms paralleled those in the 
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long-range region by employing model calculations. As an extension 
of the long-range point dipole model, Levine and McQuarrie 
represented a pair of interacting atoms by interacting metallic 
spheres. 22 The long-range asymptotic limit of the resulting 
polarizability is identical with that of the point dipole model 
presented earlier. Their short-range results indicate a sharp peak 
in the parallel polarizability at the point at which the spheres 
come in contact. The physical basis for this effect is obvious, 
since when the spheres touch, electronic charge can flow from 
one sphere to the other. They then translate this modeled result to 
mean that for real atoms overlap effects are indeed important. This 
agreement with intuition notwithstanding, their resulting value for 
B £ is positive for all cases, and greater than the point-dipole 
result by about a factor of 2. Thus one concludes that, at least 
for this problem, such classical modeling is not sufficient. 

To consider a model with internal structure, DuPre and McTague 
used the hydrogen molecule in its 3 £* excited state to simulate 
the behavior of a rare-gas diatom, 23 The justification for this 
simulation is that this first excited triplet state of H 2 displays 
the qualitative features of rare-gas diatoms in their ground states. 
That is, the H 2 potential energy curve includes a short-range 

(overlap) repulsive part with a long-range (IT 6 ) attractive part; 
this situation is characteristic of the rare-gas diatomic potential 
curves. Also, in terms of molecular orbital configurations, triplet 
hydrogen is a a , whereas singlet helium is o 2 o 2 . 

s u g u 



68 


In the DuPre and McTague calculation, a one parameter variational 
approximation to the solution of Eq„ (1.6) was used: 

Y(F) = y 0 (l + A( q ! + q 2 )) , (3.1) 

where Yq was chosen to be the Hirschfelder-Linnett wavef unction 
for the 3 J^ state of the unperturbed H 2 molecule; 52 here, qj and 
q 2 are the coordinates of electrons 1 and 2 , respectively in the 
direction of the applied field. As was shown by Hirschf elder, the 
components of the polarizability tensor of H 2 , under these conditions, 
are given by: 53 


8(<Xj> + <x 1 x 2 >) 2 , 

(3.2) 

8(<z^ + < z ^z 2 > ) 2 , 

(3.3) 


where Xj and Zj are coordinates of electron 1 in the center of 
mass system, with the z-axis lying along the internuclear axis, and 
the < > denote an expectation value with respect to y 0 . Eqs. 
(3.2) and (3,3) are a generalization to molecules of the Kirkwood 
approximation for atomic polarizabilities, 54 * 62 and were shown to 
give results of reasonable accuracy for the ground state of H 2 . The 
results for the 3 £* state of H 2 showed that the increments of a,| 
were positive and ctj were negative in the intermediate region, and 
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that in this region, the mean incremental polarizability, Aa , is 
negative* Thus, the second dielectric virial coefficient B g is 
less than that obtained from the dispersion models, and for some 
temperature range could, in fact, be negative. 

The qualitative explanation of the behavior of the incremental 
polarizabilities is based on the notion that polarizability is 
related to the size of the charge cloud. When the two atoms overlap, 
the electrons on one atom can tunnel to the other nucleus, thereby 
increasing the polarizability in the parallel direction. The same 
effect contracts the charge cloud in the perpendicular direction, 
however, and decreases this component of the polarizability tensor. 

For the case of the helium diatom, the Hirschfelder-Kirkwood 
approximation takes the form 

a !l = ^r<cr g f z 2 ] a g > +<o u |z 2 |a u > - 2<a g | z | c u > 2 ] 2 , (3.4) 

a l = 4 [ <a gf x2 l a g > + <* u !x 2 k>7 » (3.5) 

for a field-free wavefunction of the Hartree-Fock form, 

= |o a a cr | . (3.6) 

g g u u 1 

As a preliminary to our main calculation, we chose for the 
Kestner wavefunction and evaluated Eqs. (3.4) and (3.5). The 



Table III.l. Hirschfejder-Kirkwood approximation for the polarizabilities of He 2 (see Eqs. (3.4) 
and (3.5)). a,b 


R 

2<a 1 Z 2 1 a > 
g 1 1 g 

2<a 1 z 2 1 a > 

u 1 1 U 

2 <a 1 z 1 a > 
g 1 1 u 

2 <a |x 2 [a > 

g 1 1 g 

2<o |x 2 |a > 

a n 

a l 

Aa 

1.0 

0.7443 

4.7418 

1.1520 

0.4507 

1.4220 

17.100 

3.5446 

5.567 

2.0 

2.3397 

3.6952 

2.1042 

0.7797 

0.7430 

2.5832 

2,3186 

-0.089 

3.0 

4.9164 

5.7783 

3.0246 

0.8241 

0.7322 

2.3916 

2.4221 

-0.084 

4.0 

8.5505 

9.0511 

4.0042 

0.8184 

0.7546 

2.4586 

2.4743 

-0.027 

4.7 

11.6793 

11.9986 

4.7013 

0.8079 

0.7699 

2.4828 

2,4895 

-0.008 

5.2 

14.2004 

14.4096 

5.1982 

0.8021 

0.7756 

2.5240 

2,4891 

-0.005 

5.5 

15.8274 

16.0043 

5.5002 

0.7998 

0.7796 

2.4948 

2.4945 

-0.001 

6.0 

18.7322 

18.8489 

6.0003 

0.7949 

0.7841 

2.4885 

2.4932 

-0.004 

GO 

<a[z2|a> 

= 0.39494 




2.4956 

2.4956 



a All entries in this table are in atomic units. 
b The ipQ (Eq. (3.6)) is that of reference (41). 
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results are shown in Table 3 „ 1 . It is seen that this approximation 
gives negative values for Aa in the region of the van der Waals well. 

To proceed beyond these model calculations, several workers have 

applied Hartree-Fock theory. 24, 26 ,27 The first Hartree-Fock 

computation of the polarizability of the helium diatom was carried 

out by Urn, Linder, and Kromhout . 24 Their calculation used a Gaussian 

basis set and was unable to produce a negative value for the 

dielectric virial coefficient, B . However, their calculated 

e 

asymptotic results also differed from the Hartree-Fock limit for the 
separated atoms by about 25 %. Thus, their calculation is of 
questionable accuracy. 

Buckingham and Watts have also used coupled Hartree-Fock theory 
to compute an a(R) curve for He 2 using a Gaussian basis set . 26 
They obtain negative values for the incremental mean polarizability 
in the region of strong overlap. These negative contributions, 
however, occur at too short a separation to produce a negative value 
for B . Their results are plagued by instabilities in the values 
of a(R) with respect to basis set variation. To this extent, the 
authors themselves point out that their results should be regarded 
as tentative, and no claim for attainment of the Hartree-Fock limit 
is made. 

Recently, O'Brien et a},, have also published results of a 
coupled Hartree-Fock calculation of the polarizability of He 2 . 27 
Their calculation was similar in motivation and method to ours, so 



that a discussion of their results will fce deferred to the final 
section, at which time a comparison of results will be presented. 
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B. Coupled Hartree-Fock Perturbation Theory 

Langhoff, Karp! us, and Hurst have reviewed several of the forms 
of perturbed Hartree-Fock theory currently in use . 56 The method 
which was used in our work is known as the coupled Hartree-Fock 
perturbation method (method a of reference (56)). 

The Born-Oppenheimer Hamiltonian, in atomic units, for a pair 
of helium atoms relative to an origin at. the center of mass, in a 
constant, static electric field F is given by: 


4 4 

■hi vj- i (- 

1=1 1*1 i A 


2 


+ ---) 
r iB 



(3.6) 


where R ft and represent the position vectors of the nuclei. 
Invoking the Hartree-Focr approximation fcr this case involves the 
assumption that the eigenstates of H can be represented by a 
single determinant of orthoncrmal spin orbitals: 


* = dettxj.0) xz(2) X 3 ( 3 ) x 4 (4)] , (3.7) 

where the spin orbitals x are to be chosen by minimizing the energy 
functional, Eq = (1,22). The derivation and use of the resulting 
one-electron (Fock) equations is well known for the unperturbed 
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problem* 39 To extend the formalism to the problem at hand, we need 
only note that since the perturbation is a symmetric sum of 
one-electron operators, it will carry through the usual Hartree-Fock 
derivation in the same way that the nuclear attraction terms do in 
the unperturbed case* One then obtains the perturbed canonical 
Hartree-Fock equations: 


(h(l) + r v F - e.h^l) = 0 , 1 = 1,2 , (3.8) 


where h(l) is the one-electron Fock operator: 


h(l) - 


• l /iA - --- jr-+ l <x i kil (1 

1A r lB j-1 2 


P 12)|Xj > 2 > 


(3.9) 


and e ± is the perturbed orbital energy. The subscript "2" on the 
fourth term means that one integrates the indicated expression over 
the coordinates of electron 2 only. For helium, Eq. (3,8) represents 
a pair of coupled integrodifferential equations for the orbitals 
Xi and X 2 • The solution of these equations can be approached 
either directly or by using perturbation methods. 

With the direct method, one leaves the electric field explicitly 
in the Fock equations. This means that the perturbed Fock equations 
differ from the usual unperturbed Fock equations only in the electric 
field term. Thus, a standard Hartree-Fock program can be used for 



74 


the field-dependent problem with a relatively small number of 
modifications. This method was first used by Cohen and Roothaan 57 * 58 
in atomic polarizability calculations, and is known as the finite- 
field method; it is the one used in our work. 

In using the perturbation method to solve Eq. (3.8), one 
proceeds as in Chapter I, and introduces the expansions: 


, = JO) + p (1) 

r T • • - 9 


(3.10) 


and 


e i = e i 0) + F e l 1) + ••• . (3.11) 

into Eq. (3.8), and, using the usual linear independence arguments, 
one gets, to zero and first order respectively : 

D>°0) - e< 0) ] X { 0) = 0 , (3.12) 


and 


[h°0> - - [r lq + c<«]x' 0) (l) 

+ l [<X. a) | r Il (' - P 12>lx5 0> >2 (3.13) 

j=l J J 

+ < xj 0) ir-i (1 - P 12 ) Ixj 1) > 2 ]xj 0) * o, 
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where 



(3.14) 

and r lq refers to the component of in the direction of F . 
Eq. (3.13) represents a set of coupled integrodifferential equations 
for the first order contributions to the perturbed orbitals and 
energies. 

To obtain the polarizability using these two methods, we refer 
back to Eqs. (1.5) and (1.24). The field-dependent dipole moment 
for He 2 in the Hartree-Fock approximation is given by: 

<y> = <v(F)|£|?(F)> , (3.15) 

where v(F) is defined in Eq. (3.7). Denoting the orthonormal 
closed shell ot-spin orbitals for the field-dependent problem by 
Xi and X 2 » we see that Eq. (3.15) can be written as: 

K v> = 2< xilyilxi > + 2 <X 2 lyi I X 2 > • (3.16) 

This expression is used directly in the finite-field method by 
extrapolating the ratio y/F to zero field. For use with the 
perturbation method, Eq. (3.16) can be written: 
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<y> = 2<xi°^ + F xi 1 '* 4 - I y 1 1 + F 4 • ..> 

+ 2<X2 0) 4 F X2 15 + ••.|jjilx2°^ 4 F xP 4 ...> (3 o 1 7 ) 

= 4F[<xi 0) Ui|xP> + <x| 0) I 'Ml I x| 1);> Il 4 0(F 2 ) . 


Here, we have used fact that the dipole moment of the free He 2 diatom 
is zero. Differentiating Eq. (3,17) with respect to F , and 
evaluating the result at zero field yields the following 
expression for the polarizability: 


a 


qq 


2 


4 l 
1=1 





(3,18) 


In both the finite-field and perturbation methods, exact 
solution of the Hartree-Fock equations is not in general possible, 
so that one makes use of basis set expansions to obtain analytical 
solutions. That is, in the finite-field method, one uses the 
Hartree-Fock-Roothaan procedure by introducing a basis set expansion 
for the field-dependent orbitals x ± . This converts Eq, (3,8) to a 
set of coupled matrix equations for the expansion coefficients. In 
the perturbation method, one converts to a matrix scheme by 
introducing a basis set expansion for the first order orbitals, 
xP . The resulting coupled matrix equations are solved iteratively 
in either case. It is readily seen that for small field strengths, 
the two methods are equivalent. 
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In both the finite-field and perturbation methods, one begins 
with a zero-order Hartree-Fock wavefunction which, except for certain 
atomic cases, is itself obtained by the Hartree-Fock-Roothaan 
procedure. The basis sets used in such calculations usually consist 
of either Slater type functions (STF) or Gaussian type functions 
(GTF). The STF sets are generally more efficient than the GTF, 
with the size of the latter often being three or more times larger 
than the former for the same accuracy. 

The symmetry of the zero order basis functions is dictated by 
the molecular orbital correlation diagram for the particular 
molecule of interest. 59 The set of functions which one introduces 
to reflect the distortion of the zero order wavefunctions by the 
electric field is known as the polarization set. For polarizability 
calculations, the symmetry of these functions is chosen to be such 
that they couple to the zero order molecular orbitals via the dipole 
moment operator. Sitter and Hurst have this discussed this problem 
for calculations of atomic hyperpolarizabi 1 ities . 60 

In view of the expansion nature of the solutions of Eq. (3.8), 
it is important to ask how the result obtained compares to that 
which would be obtained if the Hartree-Fock equations could be solved 
exactly. This latter limit is referred to as the Hartree-Fock limit, 
and it corresponds to the use of an infinitely large basis set. The 
difference between it and the analytical result is known as the basis 
set error, and is exactly known only for atomic calculations. The 
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"art" of carrying out such computations, therefore, involves the 
minimization of the size cf the basis set used, while, at the same 
time, getting as close as possible to the Hartree-Fock limit result- 
For field-free diatomi molecules, this problem is considered to be 
solved, and Gilbert and Wahl have illustrated it in detail for the 
He 2 , Ne 2 , and Ar 2 rare-gas diatoms- 61 Additional references to this 
problem and its solution for field-free calculations can be found in 
Chapters I and III of reference (63). 

The method which we used fo*“ the selection of the perturbed 
basis set will be discussed in the next section. 

From a formal point of view, the finite-field procedure Is the 
same whether one is discussing atoms or molecules. From a 
computational standpoint, however, this is no longer true because 
of the additional amount of work that one must do in computing and 
processing the requisite two-center, two-electron integrals. For 
this reason it is desireable to recast the above procedure to take 
advantage of optimized computer codes which currently exist. 

The program which we used in our calculations is a modification 
of the BISON computer program of A. C- Wahl, P. J. Bertoncini, R. Land, 
and K. Kaiser of Argonne National Laboratory. 6 ^ The modifications are 
those required for field-dependent calculations, and are discussed in 
more detail in Appendix C. 

The basic change needed involved the introduction of the 
appropriate products of the dipole moment matrix elements and the 
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electric field strength in the one-electron contributions to the Fock 
matrix. Since the field variation was done in the SCF part of the 
program, unit field strengths were used in the one-electron integral 
calculations; this means that the integrals need to be computed only 
once per each R value considered. This was especially important 
in that convergence problems which appeared in the calculations of 
ot|| (R) , made it necessary to turn the field on slowly. 

For basis set variation purposes, auxiliary programs were 
written to remove those one- and two-electron integrals which 
involved a given basis function from a "master" integral tape. This 
allowed step-by-step basis set reduction to be performed using the 
Integrals corresponding to a single, large set. One, of course, 
must be careful to guarantee that the cost of computing the integrals 
for the large set is less than the sum of the costs of the individual, 
smaller sets. 

Once the appropriately modified Hartree-Fock program was working, 
the problem remained of delineating a basis set selection procedure 
which would allow us to approach the Hartree-Fock limit economically. 
The procedure which we used is discussed in the next section. 

C. Basis Set Selection Procedure 

In order to compute the Hartree-Fock limit for the polarizability, 
one introduces basis functions which reflect the electric field 
polarization of a zero-field Hartree-Fock limit wavefunction and 
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carries out the Hartree-Fock-Roothaan prescription with this augmented 
basis set. This procedure for the polarizability calculation can be 
capsulized as follows: 

1. Obtain a field-free Hartree-Fock limit wavefunction using 
the Hartree-Fock-Roothaan method. 

2. Choose polarization functions for the above. 

3. Redo the Hartree-Fock-Roothaan calculation with the zero- 
order plus polarization basis set in the presence of an 
electric field F . Compute the dipole moment y of the 
resulting field-dependent wavefunction and extrapolate the 
ratio y/F to zero field to get a . 

We now assume that step 1 has been carried out. The choice of 
a polarization set for a given field-free basis is a playoff between 
two competing requirements: 


1. The set must be flexible enough to adequately represent the 
distortion of the field-free orbitals by the electric field. 

2. The size should be minimized, within the above constraint, 
so as to minimize the number of integrals which must be 
computed . 

The prescription which we used for the selection of polarization 
functions will now be described. 


As a first step, the polarization functions are chosen to 
polarize individual field-free basis functions. For a given 
normalized Slater orbital, we define the orbital polarizability 

by the maximum value of the second-order variational functional: 
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a($ (l) ) = -2<$ (1) |(H 0 - E 0 ) + 4«f> (1) |zj<t> (0) > , 

(3.19) 

where is a trial polarization function, z is the polar axis 

of coordinates, and H 0 is an unperturbed Hamiltonian which defines 


H 0 = -^V 2 - 

tnr -1 - [jiU+ 1) - n(n-l)]/2r 2 , 

(3.20) 


E 0 = -S 2 /2 > 

(3.21) 

AO) - 
*0 

Nr n-1 e -Sr Y te (e,4) . 

(3.22) 


We expand as a linear combination of STO's and maximize 

ct($^) with respect to the linear coefficients and the orbital 
exponents. A general nodeless STO is the ground state of the 
corresponding H 0 , so that for a sufficiently flexible choice of 
, the maximization procedure yields the exact polarizability 
of the ground state of H 0 . Following the optimization of , 

we generate a more flexible basis of polarization functions by 
adding the STO's which represent , to the finite-field basis 

individually; that is, we use the variational procedure to choose 
optimum orbital exponents and principle quantum numbers for the 
finite-field basis. 
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If the procedure described in the preceding paragraph is 
followed to polarize an extended field-free basis, a linearly 
dependent basis set is usually obtained. Even if the linear 
dependencies are removed, we have found that the basis is much 
larger than is necessary to give converged results for the molecular 
polarizability. We then consider the following modifications to 
reduce the size of the polarization basis: 

a. Primary attention is paid to polarizing the STO's which 

fall into the following categories: those which are dominant in the 

field-free orbitals; those which are the most polarizable; and those 
which constitute a leading contribution to the tail of the orbitals. 

b. The less dominant basis functions are polarized collectively. 

Let and denote two STO's which have expansion 

coefficients Cp and c^ in a field-free occupied orbital, and let 

and <t>^ denote their polarization functions. Then we 
construct a new polarization function 6 , which also is represented 
by an expansion in STO's, in a weighted least-squares sense. We 
minimize 

a = | |e(r) - £ (1) (r)[ 2 dx , (3.23) 

where 
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4 


a 


p 



c <f>^ + a c , 

p p q q q 

(3.24) 


(3.25) 


The optimum expansion functions of e give a smaller set of STO's 

to add to the polarization basis than the two functions d>^ and 

p 

taken individually. 

Other factors affecting the choice of basis functions depend on 
the case under consideration and are discussed below. 


D. Test Calculations 

Before computing the polarizability of He 2 , we obtained the 
polarizabilities for the united and separated atoms, Be and He, 
respectively, in order to test our procedure for basis set selection; 
the idea being to obtain the coupled Hartree-Fock limit for the 
polarizability of these atomic limits. The results for the atomic 
tests are given in Table III. 2. These results agree quite well with 
the accepted coupled Hartree-Fock values (a(He) = 1.322 a^ ; 
a(Be) = 45.5 a^). 57 * 58 ’ 65 

The atomic basis set selection followed closely the procedure 
outlined in the previous section. The field-free basis for Be is 
taken from Clementi 66 while that for He is taken from Kestner. 41 The 
3p orbital exponent in the He atom calculation was selected by 
variationally solving Eq. (3.19), for each 2s basis orbital, using a 
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two-term basis set (2p,3p) for ; the combined result was then 

fitted to a single 3p via the least-squares procedure (see Eq. (3.23)). 
For the Be atom, the selection of candidates for collective 
polarization is not clear-cut, because the second occupied orbital 
has a dominant contribution from the 2s orbitals, while the first 
has a dominant contribution from the Is orbitals. Thus, these 
orbitals were polarized separately. 

In addition to these atomic limits, a series of test calculations 
were performed on the hydrogen molecule, at an internuclear separation 
of 1,402 a 0 . These tests were made both to check the modified BISON 
program and to examine features of the finite-field method and our 
basis set selection procedure peculiar to a diatomic example. 

The program checks were of two types: the first involved a 

comparison with a coupled Hartree-Fock calculation which was carried 
out by hand; and the second involved the duplication of results 
obtained by other workers who used the iterative method. 

The former check is possible without iteration because H 2 has 
an electronic configuration which can be represented by la* t s0 that 
there is only one perturbed orbital, la (1) , hence only one equation 
of the form of Eq. (3.13) to solve. This means that expansion of 
la^ ^ in a Slater basis set leads to a single matrix equation of 
the form: 


b , 


Ac 


(3.26) 
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which is solved by inverting the matrix A . The second method 
involved the duplication of the coupled Hartree-Fock H 2 polarizability 
calculations by Caves and Karplus 67 and Gutschick and McKoy. 68 

For the first test, Eq, (3.13) for the field perpendicular to 
the internuclear axis can be written: 


C-^ v i - - -pjj + <0 o I I a o > 2 - 

+ Cxi + 2<a 0 |r~2|a l > 2 - e 0 ]cfg = 0 • 


(3.27) 


Here, ej is given by (aoailffQOo) . If we now introduce the 
expansion. 


01 “ l c K * K . 


K 


(3,28) 


into Eq. (3„13) we obtain Eq. (3.25), with. 


XK 


+ 2 ( ctq 4> l } CTo<t> K ) , 


(3.29) 


and 


b L = - <a olxU L > • 


(3.30) 
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The required integrals were computed using the two-center integral 
program of H. F. Schaefer I I I „ The results of this check for 
are given in Table III. 3. The minor differences reflect the (arbitrary) 
convergence tolerance of 10“ 5 for the coefficient vectors in the SCF 
part of BISON. This check lends confidence to the correctness of our 
working version of BISON. 

To further check the program, we used the basis set of Caves 
and Karplus, and computed a y and a, for H 2 at R = 1.40 a 0 . 

The results are given in Table IIL4. Again, the slight discrepancies 
result from differences in numerical parameters used in the 
respective calculations. These results serve to re-enforce the above 
conclusions regarding the accuracy of our working version of BISON. 

A further concern which one has in using the finite-field 
method is that hyperpolarizabil ity terms can contribute to the 
computed value of the dipole moment. If these additional effects are 
not correctly interpreted and accounted for, erroneous results for 
the polarizability will be obtained. 

For a calculation in which the polarizability is the only 
quantity of interest, the hyperpolarizability terms can be removed 
by using a sufficiently small field strength. However, the field 
strength must be kept large enough so that the computed dipole 
moment has numerical significance. The particular value of the 
field strength used in these calculations was .001 au . This value 
was chosen by test calculations on both H 2 and He 2 . 
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Table III. 2. Basis sets and polarizabilities for atomic helium and 
beryllium.. 


Atom 

Field-Free Basis 9 

Polarization Basis 3 

Polarizability 13 


1 s ( 1 .450) 

2p( 1.450) 


He 

2s(l .641 ) 

3p(l .450) 

1.322 a 3 
0 


2s(l .723) 

3p(l .655) 



1 s (3 . 3370) 

2p,3p(3.3370) 



1 s ( 5 . 5063 ) 

2p,3p(5.5063) 

* 

Be 



45.574 a! 


2s(0.6040) 

2p,3p(0.6040) 

0 


2s(l .0118) 

2p,3p(l .0118) 



cl 

Slater- type basis. Orbital exponents given in parentheses. 
b Here a 0 is the Bohr radius. 



Table III. 3. Comparison between the first order perturbed orbital for H 2 at R = 1.40 a 0 computed 
using Eq. (3.27) and that obtained from BISON. 3 


An" 

A 12 

A 22 

bi 

b 2 

F 

Cl 

C 2 

Eq. (3.27) 2.4111 

BISON 

2.4658 

3.5624 

-1.1667 

-1.0919 

0.025 

0.025 

-0.01458 

-0.01457 

0.00243 

0.00243 

3 In the basis set definitions, the 

subscripts A and B 

refer to 

the respective 

centers: 



cj q — -3.62184>i ” .1 7759<}i2 ” .O 28 I 743 

♦ 1 “ XiqqO-06)^ + xiooO-06) b 
<f>2 = XiooO ’43)^ + XlooO‘43)g 
$3 = X21o(1-82) a + X21o(^*82)g 
a l = C I $4 + C'l 45 

H = Cx2ii(.935) A + X2I - l (.935) A + X2n(*935)g + X 2 1 — 1 C »935)g] 

<J>5 ~ Cx2ii(l»509)^ + X 2 i-i(l ”509)^ + X2ii(^*^09)g + X 21- 1 ( 1 ° 509)gJ . 

The x's are defined by Eq. (3.22). 
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Table II I. 4. Coupled Hartree-Fock polarizabilities and basis sets 
for H 2 at R = 1.40 . 


Field-Free Basis' 5 

Polarization Basis 3 

Polarizability 

ls(l .06) 


6.478 a 3 

a.. ls(l .43) 

2p (0.945) 

U 

II 

2p ff (1.82) 

0 

(6.478 a 0 3 ) b 



1 s (1 .06) 

2p (0.935) 

4.601 a 3 

a l 

1 s (1 .43) 

2p o (1.82) 

IT 

2p (1.509) 

7T 

U 

(4.602 a 3 ) b 


a Numbers in parentheses are the orbital exponents used, with a 


Slater-type basis centered on the respective nuclei, 
b 


Reference (67). 
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Table III. 5. 

Coupled Hartree-Fock ratios at 
R = 1.402 a 0 . a 

<y x >/F , for H 2 at 

F 

"x< F >I b 

"*< F >n c 

0.005 

4.290 

4,290 

0.008 

4.289 

4.290 

0.010 

4.289 

4,290 

0.025 

4.281 

4.292 

0.050 

4.253 

4.294 


Units of field strength are atomic units, units of y /F are 

X 


1 refers to the basis set: ls(1.197), 2p ir (1.050), 3p (1.197) on 

each center. 

c II refers to the above basis set (b) augmented by the function 
3d A (1.197) on each center. 
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The hyperpolarizability of a diatomic molecule will depend on 
basis functions of 2 , tt , and A symmetry. Thus if one hopes to 
extract hyperpolarizability information from a finite-field 
calculation (this capability makes the finite-field method 
attractive), basis functions of all the appropriate symmetries for a 
given direction of the applied field must be included. Failure to 
do so can result in the wrong sign for the hyperpolarizability, as 
evidenced by the H 2 data of Table III. 5. 

In addition to the tests described above, we repeated the test 
calculations of Gutschick and McKoy on H 2 , and verified their results. 

E. Helium Diatom Calculations 

For our calculation of the polarizability of the helium diatom, 
the field-free wavefunction of Kestner was chosen. 41 This 
wavefunction is defined in terms of a Slater orbital basis set, with 
orbital exponents optimized at each internuclear separation. The 
basis set error in Kestner's calculation, as inferred from the 
extremely accurate calculation of Gilbert and Wahl, is about 1 part 
in 10 7 for the energy. A reduction in the size of this zero-order 
set can be justified only by numerical experimentation. 

For the purpose of calculating statistical averages, one is 
interested in the polarizability for those internuclear separations 
for which the statistical pair distribution function is non-zero. 

As will be seen later, this implies that the smallest internuclear 
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separation for which we need to compute a is about 4a o . Thus, 
for this range of separations, it is reasonable to attempt to achieve 
a reduction of the size of the basis set by experimentation at 
R = 4a 0 , and assume that the results are valid for larger separations. 

In the H 2 test calculations, it was observed that it was 
unnecessary to include basis functions which polarize the 2p 

0 

zero-order basis function. This is presumably due to two effects: 

the 2 p ct does not easily distort in a direction perpendicular to 

the orbital axis, and the coefficient of the 2p in the zero-order 

0 

occupied orbitals is small relative to those of the s-type basis 
functions. For the case of He 2 , the coefficient of the p-type basis 
functions is less than .5% of the coefficients of the s-type functions, 
for internuclear separations greater than or equal to 4a 0 . Thus, it 
was assumed at the outset that for the calculations of , the 2p 
need not be polarized. 

The zero-order basis set for R =-* 4a 0 , is given in Table III. 6. 
From this table, it is seen that for a calculation of 
a reasonable initial polarization set would consist of the following 
functions on each center: 

2p ± (1.458) , 3p ± (l .458) , 2p ± (1.175) , 3p + (1.644) . 

The last two polarization functions represent collective polarization 
of the two 2s zero-order basis functions. Since the zero-order set 
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Table III. 6. Field-free basis set for He 2 at R = 4.0 ag , 

Basis 3 Field-Free Basis 

ls(l .458) 

2s (2. 631) 

1 

2s (1. 723) 

2p o (1.748) 

3 Slater-type functions centered on nuclei A and B . Orbital 
exponents given in parentheses. 
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consists of eight functions, this yields a total of 24 basis functions, 
or 45,150 two-center, two-electron integrals which need to be computed., 
The systematic reduction in the basis set size for R = 4a 0 is given 
in Table III. 7. The most noticeable conclusion warranted by these 
results is that in the calculations of a . , the 2p basis function 
does not have to be included. This result coupled with the collective 
polarization of the two 2s functions by a single 3p yields a basis 
set of 18 functions, or 14,706 two-center, two-electron integrals. 

This corresponds to a 65% decrease in the cost of computing each point 
on the a^{R) curve (with calculations done in double precision on 
the UN1VAC 1108). 

For the parallel direction, a set of 18 basis functions was used 
for size-consistency with the calculations for the perpendicular 
direction. The size here represents a linear independent combination 
of an essentially complete polarization set; for the purposes of this 
calculation, the linear independence criterion which was adopted was 
that the smallest eigenvalue of the overlap matrix could not be less 
than 10-6 , The specific polarization sets and the results for the 
polarizabilities are given in Table III. 8 ( 0| ) and Table III. 9 ( aj ). 

The mean incremental polarizability, & a , which is computed from the 
data in these tables is plotted in Figure I II , 1 . The cross-hatched 
area in this graph reflects our estimates of the basis set error in 
the polarizability calculations, and will be discussed in Chapter V. 

From the data of Tables 1 1 1. 8 and 1 1 1. 9, it is apparent that, 
in the overlap region, a |( first increases and then decreases as 
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the Inter-nuclear separation decreases,, This is opposed to the trend 
for the corresponding values which decrease monotonically with 
decreasing R . Hence, our results indicate that the incremental 
parallel and perpendicular polarizabilities have opposite signs for 
the range of internuclear separations which we have considered. 

This, as noted in references (23) and (37), is in accord with the 
general features of the polarizability curves which are necessary to 
obtain agreement with the dielectric constant experiments. The Kerr 
effect experiments, on the other hand, are affected by the variation 
in the anisotropy as a function of R . Our computed B(R) curve 
is given in Figure III. 2, and will be discussed later. 

In order to cast this data in a form which is useful for the 
property calculations, we fit the incremental mean polarizability, 
Act(R) and the anisotropy p(R) to the following forms: 


MR) = -12.9exp(-1.45)R + .00116a* ] 


(3.31) 


and 


f R <_ 6.5a 0 . 


P(R) = aj| - = 6ctQ R -3 - 42.5exp(-1 .645R) 


For internuclear separations greater than 6.5a 0 , we used the 
asymptotic results: 
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Table I II. 7, Basis sets and perpendicular polarizabilities for 
diatomic helium, R - 4.0 a 0 . 


Basis 3 

Field-Free Basis 

Polarization Basis 

Polarizability^ 


1 s ( 1 .458) 

2p + (l .458) 


1 

2s (2. 631) 

3p ± { 1.458) 

2.570 a„ 


2s(l .723) 

2 P± (1J75) 



2p o (1.748) 

3p ± (l .644) 



1 s ( 1 .458) 

2p ± (l .458) 


2 

2s{2.631 ) 

3p + (l .458) 

2.570 a* 


2s 0 - 723) 

2p ± (l .175) 

V 



3p ± (l .644) 



ls(l .458) 

2 P± (1 .458) 


3 

2s{2.631) 

3p + (l .458) 

2.570 aj 


2s(l .723) 

3p ± (l .644) 


Slater-type basis centered 

on nuclei A and B . 

Orbital 


exponents given in parentheses, 
k Here a 0 is the Bohr radius. 
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Table III. 8. Basis sets and parallel polarizabilities of diatomic 
helium. 


R(a 0 ) -*• 

Basis 

I 

4.0 

4.7 

5.5 

6.0 

oo 



Orbital Exponents 



Is 

1 .458 

1.454 

1 .450 

1.451 

1.450 

2s 

2.631 

2.635 

2.630 

2.638 

2.641 

2s 1 

1.723 

1.723 

1.721 

1.720 

1.720 

2p 

1.250 

1.250 

1.250 

1.250 

1.250 

2 P* 

1.458 

1.454 

1.450 

1.451 

1.450 

3 Pa 

1.458 

1.454 

1.450 

1.451 

1.450 

3p 

1.655 

1.655 

1.655 

1.655 

1.655 

3d 

a 

1 .458 

1.454 

1.450 

1 .451 

1.450 

4d 

a 

1.458 

1.454 

1.450 

1.451 

1.450 


Parallel Polarizabilities (ajj) 


(a) 

2.675 

2.685 

2.681 

2.672 

2.644 

(b) 

2.678 

— 

2.681 

2.675 

2.650 

(a) This work. 

Here 

a 0 is the Bohr radius. 




(b) Gaussian basis set results of reference (27). 
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Table III. 9. Basis sets and perpendicular polarizabilities of 
diatom helium. 


Basis 

4* 

R(ag) -»■ 

4.0 

4.7 

5.5 

6.0 

CO 




Orbi tal 

Exponents 



Is 


1.458 

1.454 

1.450 

1.450 

1.450 

2s 


2.631 

2.635 

2.630 

2.641 

2.641 

2s' 


1.723 

1.723 

1.721 

1.723 

1.723 

2p ± 


1.458 

1.454 

1.450 

1.450 

1.450 

3p ± 


1.458 

1.454 

1.450 

1.450 

1.450 

3 P± 


1.6550 

1.6550 

1.6550 

1.6550 

1.6550 



Perpendicular Polarizabilities (a q ) 


(a) 


2.570 

2.602 

2.623 

2.628 

2.644 

(b) 


2.582 

— 

2.630 

2.635 

2.650 

(a) 

This work. Here 

a 0 is the 

Bohr radius. 


(b) 

Gaussian 

basis set 

results of reference 

(27). 
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( 


Figure III.l. Mean Incremental Polarizability, Aa(R) , for He 2 , 
as a Function of R . 



R (a 0 ) 
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Figure III. 2. Anisotropy in the Polarizabil ity Tensor, p(R) , 
as a Function of R . 
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Aa(R) = 4a o R- 6 ) 


(3.33) 


and 


f R > 6.5a 0 . 


P(R) = 6al R" 3 + 12.07R- 6 J 


(3.34) 


The successful use of an exponential function to fit our 
computed mean incremental polarizabilities in the overlap region 
bears no special physical or mathematical significance. In fact, 
attempts at fitting the computed parallel and perpendicular curves 
individually, using a similar functional form, were unsuccessful, so 
that one cannot justify the use of an exponential fit to our data on 
any grounds other than convenience. Here, a 0 is the computed 
atomic polarizability of helium (1.322ao) . The use of the long-range 
extrapolation is justified in the next section. 


F. Extrapolation of the Hartree-Fock Polarization to the 
Long-Range Region 

Our calculations do not extend beyond an internuclear separation 
of 6a 0 because the increment in the polarizability compared to the 
atomic limit becomes zero to the accuracy of our calculation. 

However, the asymptotic atomic limit of the coupled Hartree-Fock 
polarizability can be derived without explicit calculation. 
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It is well established computationally that the Hartree-Fock 
potential energy curves for interacting closed shell atoms are 
strictly repulsive and do not account for van der Waals (R -6 ) 
binding. Qualitatively this occurs because the effective potential 
for the localized orbitals on one atom is obtained by spatial 
averaging of the instantaneous multi pole moments of the other atom. 
Since the atoms have no permanent multipole moments, the interactions 
average to zero, leaving only the repulsive overlap forces as 
nonvanishing. 

In considering the incremental polarizability, a similar 
situation occurs with the uncoupled Hartree-Fock method. 56 ’ 69 Here 
one calculates the polarizability of the field-free orbitals 
separately, without allowing the remaining occupied orbitals to 
simultaneously polarize. Hence, at long range where overlap effects 
can be neglected, the polarizability of each localized atomic orbital 
in the diatom is equal to its atomic value, and there are no 
multipole contributions from the other spherical atom. The 
asymptotic limit of the uncoupled theory is then, for either 
component oc|| or , 

a(uncoupled) - 2a 0 + overlap terms , 


where a 0 is the uncoupled Hartree-Fock polarizability of the atom. 



105 


A more interesting situation occurs in a coupled theory, 58 ’ 69 
such as we have used. Here the polarizability of each orbital is 
computed with the orbital in the presence of the remaining polarized 
orbitals. Thus, at large internuclear separations, a localized 
orbital on one atom feels not only the external field, but also the 
field due to the induced dipoles on the other atom. This picture 
is reminiscent of the classical point-dipole model of the incremental 
polarizability. In fact, we can show that the long-range limit of 
the coupled Hartree-Fock theory gives precisely the point-dipole 
model results; namely, 

a,, - 2a 0 + 4a § R~ 3 + 8a£ R“ 6 + ... , (3.35) 

~ 2a o - 2a. q R ~ 3 + 2<xg R~ 6 + . . . , (3.36) 


where a 0 is the coupled Hartree-Fock polarizability of an isolated 
helium atom. (We restrict attention to helium, although the results 
of this section are general for interacting closed shell atoms.) 

To prove Eqs, (3.35) and (3.36) we consider two helium atoms at 
large separations R . The Hartree-Fock energy of the diatom in the 
presence of an external field F may be written 


E - 


2h + J 

a aa 


2u *F + 2h, + J 

-a ~ b bb 


2u b *F' + 4R- 3 (y -y 


3y y, ) 
az bz / 




(3.37) 
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where we have used a localized representation of the wavefunction, 
ip - a(l )a(2)b(3)c(4) , and the orb<ta- a is centered on nucleus 
A . We have retained only the leading (R" 3 ) multipole interaction. 
The various terms in Eq. (3.37) are expectation values of the 
corresponding operators; e.g. 

y < a I y | a / a | a> » (3.38) 

^3. '3 


and 


J _ * <a|J |a i a|a> , (3.39) 

3LSL d 

A 

where J is the Coulomb operates 

J =■ | a 2 (2)r-^ c! , ' i a ? ( ? ) 6-2 - (3.40) 

The Fock equation for the orbita’ a Is obtained by varying 
the energy E with respect to a » while holding the orbital b 
fixed. Requiring the energy to be stationary results in the equation 


[h + 
a 


J - y »F + 2 R" 3 (l 


y *r -r £K »\y .y - 3 U y 
a ~a ~ tb ta b z' a z 


) - f]a - 0 


(3.41) 


The proof of the point-dipole model now rests on showing that a self- 
consistent solution of Eq. (3.41) has the form of a scaled atomic 
orbital , 
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a “ a 0 (cF) , (3.42) 

where ? is a constant to be determined. Here the orbital a 0 
satisfies Eq. (3.41) for R =■ » . By symmetry a similar solution 
holds for the orbital b . 

Let us consider two geometries. First suppose F is 
perpendicular to the internuclear (z) axis; then Eq. (3.42) for b 
gives 

y b = <b|p b Jb>/<b|b> = p(?F)?? . (3.43) 

The leading term in p is one-half the atomic polarizability a 0 , 
while higher-order terms in F give hyperpolarizabilities. 
Substituting Eq. (3.43) into Eq, (3.41) we see that the effect of 
the multipole term is simply to scale the external field. Thus the 
condition that Eq. (3.42) be a solution to Eq. (3.41) is 

4 * [1 + 2p(cF)R -3 ] -1 . (3.44) 

With this result for t, , the dipole moment of the diatom is 

2y a + 2u b = 4p(l + 2p R" 3 ) -1 F . 
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The diatomic polarizability is the leading term of the expansion of 
the dipole moment in powers of F ; thus, 

= 2a 0 (l + a 0 R“ 3 )"l . (3.45) 

This is the point-dipole result. 

Similarly, an external field parallel to the internuclear axis 

gives 


“a = 2o 0 (l - 2a 0 R" 3 )" 1 , (3.46) 

which is also the point-dipole result. 

If higher multipoles are retained in Eq. (3.37), more complicated 
expressions involving quadrupole polarizabilities, etc. are obtained. 
This also occurs in the classical point-dipole model if field- 
gradient effects are considered. Such terms do not affect the 
polarizability through order R“ 6 , however. 

In our computations using the finite field method, we did not 
achieve sufficient numerical accuracy in the long-range region 
(R > 6a 0 ) to observe the asymptotic limits. Nevertheless, in using 
the computed polarizability functions, it is legitimate to 
extrapolate to large internuclear separations via the point-dipole 
model . 


CHAPTER IV 


STATISTICAL METHODS 


In Chapter I it was noted that, in addition to the polarizability 
curves, ot u ( R ) and a^(R) , the computation of the properties in 
which we are interested require a knowledge of the function g(R) 
for gaseous helium. This function is the density independent part 
of the first term in a cluster development of the pair distribution 
function; that is, it is the pair distribution function for a dilute 
gas. For high temperatures, this simply involves a knowledge of the 
interaction potential U(R) for a pair of atoms. In the low 
temperature regime, however, things are not quite so easy, and one 
needs to compute the so-called Slater sum, W 2 (q) , where q 
represents a complete set of generalized coordinates. 

The differences between the classical and quantum mechanical 
cases arise because the state of a quantal system is specified by a 
wavefunction, rather than the exact values of the coordinates and 
momenta, as it is in classical mechanics. In addition, it follows 
that spin effects give rise to terms in the quantum mechanical g(R) 
which have no counterpart in the classical expression. In fact, the 
quantum mechanical g(R) can be written as a sum of two terms: 
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no 


g(R) g ( R )direct 1 9 ^exchange * 

where the first term is the quantal g(R) for a Boltzmann gas, and 
reduces to the classical Boltzmann factor in the limit of high 
temperature. The second term represents spin effects, and is 
associated with the Bose-Einstein (+) or Fermi-Dirac (-) 
character of the gas. A discussion of this latter point has been 
given by Yang and Lee, 70 and several workers have illustrated it with 
mode" calculations. Larsen has computed the two terms of Eq. (4.1) 
for a gas of hard spheres; 71 Larsen, Witte, and Kilpatrick have 
obtained results for 4 He at low temperatures using a Lennard-Jones 
12-6 pair-potential; 72 while Poll and Miller also used a 12-6 
potential to compute g^ and 9 exc ^ an g e for H 2 gas for 
temperatures in the range of 2°K to 80°K , 73 

The method used in this work to compute g(R) is that of Klemm 
and Storer. 74 > 75 For high temperatures, one can approximate the 
statistical density matrix P (R,R' ;b) (e = 1/kT) , by expanding it 
in powers of @ , and neglecting terms which are greater than second 
order in e . The density matrix, p , is then given by; 

p = Pi p ( 0) Pi , (4.2) 

where the operators pj and satisfy the differential 

equations 



m 


M h ip i 


9p i 

a? - ’ 


(4.3) 


and 


H o p<0) * - 2£ ?r > (4.4) 

respectively, with boundary conditions 

lim P C0) (R,R';e) = 5(R - R') = 11m Pl (R,R';e) . (4.5) 

e + o ~ ~ ~~ e->o 

Here, Hi and H 0 refer to a partitioning of the Hamiltonian of 
relative motion, H r , for a pair of particles. It is generally taken 
to be the case that H : contains the interatomic potential term U(R) 

of H r . 

Once this high temperature limit has been obtained, one can 
proceed to lower temperatures in a recursive manner, using the 
following formula due to Storer: 76 * 77 

p(R,R';2e) = | dR 2 p(R,R 2 ;6)p(R 2 ,R';3) dR 2 , (4.6) 

which is a matrix form of the operator relation 


e-2BH = e' 6H e- BH 
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If the assumption that the interatomic potential U(R) is 
spherically symmetric is made, the eigenfunctions of H r can be 
expanded in terms of spherical harmonics, so that the density matrix 
is expanded in the following manner: 


p(?.R';e) = l |^p a (R,R';b) Pjcosw) . (4.7) 

£ 0 

where P £ is the Legendre polynomial of order z , w is the angle 
between R and R’ , and p^(R,R';&) is defined by the 
differential equation. 


2P 

H2 ae 


3R2 


[ M r i + #v(R)]p 2 - 


■«£ P * * 


(4.8) 


where y is the reduced mass. The boundary conditions are: 


lim P (R,R';e) = S(R-R') for R,R' > 0 , (4.9) 

3*0 


and 


p a (R,R'; 0) = 0 for R,R' = 0 . (4.10) 

These equations can be transformed to be compatible with the form of 
Eqs. (4.3) and (4.4), by defining, 



* 


(4.11) 


H 


0 

i 



+ iL^il 

R2 


and 


= U(R) , 


(4.12) 


so that the analogues of Eqs. (4.3) and (4,4) are respectively: 


and 


fcU(R)Pi = - 


92 p 1 0) &U+1) _ (Q) _ 3p | 0) 

3R2 " R2 p j> 3g 


(4.13) 


(4.14) 


The solutions of this latter equation are the t-th partial wave 
contributions to the free particle density matrix. 


p< 0) (R,R'; B ) - 4*RR' exp[-(R2 + R'*)/(^)]y§^> , 

(4.15) 

where i is the Modified Spherical Bessel function of the first 
kind, defined by: 70 
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Thus, we can make the initial high temperature approximation: 


p £ (R,R' 


;e) = 


( ^ )U(R,) p, (0) (R,R' ;B )e-f^) u ( R ) , (4.17) 


and extend this lower temperatures using Eq. (4.6), which for the 
partial wave expansion becomes: 


p a (R»R' ;2e) 



P^{R»R2>6)p^(R2»R' >&) 


dR 2 


(4.18) 


In practice, one evaluates the various functions over a finite grid, 
so that Eq. (4.18) becomes an n x n matrix. Thus one can use a 
Simpson's rule quadrature to evaluate the right of Eq. (4.18) to 
reduce it to a square of the n x n density matrix. One then 
proceeds to lower temperatures by performing standard matrix 
multipl ication. 

The expression for g d1r and 9 exc j iange for a given value of 
6 , in terms of the partial wave expansion are: 


g d 1 r (R) 


v y ' 



(4.19) 


and 


9 exchange^ = (2S+1) -2 p £ (R,R;6) . 


(4.20) 
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A detailed account of this method and a listing of the computer 
program which is used to carry it out, can be found in reference (75). 

For the helium isotopes 3 He and 4 He (the masses used in this 
work are respectively 3.0160 and 4.0026 amu ) low temperature 
calculations of g(R) were made at T - 4°K . To get down to this 
temperature, a starting temperature of 1024°K was used, so that 
eight iterations were required. The adequacy of this starting 
temperature was established by comparison of our 4 He results at 2°K 
(9 iterations) with those of reference (72); agreement to at least 
three places was obtained. In this work, we used a mesh size of 
.025 (in units of o - 2.556 A ) at higher temperatures, and .10 
at lower temperatures for the Simpson's rule integration of Eq. (4.18). 
By numerical experimentation, we found that twenty partial waves were 
sufficient for three decimal place accuracy. In this work, the 
MDD-2 pair potential of Bruch and McGee was used. 12 This potential 
combines a Morse potential for separations less that 6.96a 0 , with 
a dispersion dipole plus quadrupole result for separations greater 
than or equal this core. The potential parameters for helium along 
with the resulting asymptotic forms of g (R) are given in Table IV. 1. 

For future reference, a plot of the values of g(R) for 3 He 
and 4 He at T = 4°K is given in Figure IV. 1, and the g(R) values 
are listed in Table IV. 2. 
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Table IV. 1. MDD-2 pair-potential and long-range g(R) for 
helium. a ’ bjC 


MDD-2 Pair-Potential: 

1 w/n\ . 36.4778 15.0935 

F VCR) p p 


1 V(R) = 10.75[e^ 2c ^" x ^ - 2e^ 1_x ^] d 


R > 1,44085 o- units 


R < 1.44085 a-units 


Long-Range g (R) Expressions: 


g(R) 

= 1 - 

V(R) 

■ft 2 v e p 

> 3,1 o- 

kT 

6mk 2 T z V R 

3 He,T = 4°K 

g(R) 

= 1 

l 9.11945 L 31.8405 . 
+ ~R5 + 

21.6784 
1 R t & 

4 He,T = 4°K 

g(R) 

<= 1 

. 9.11945 24.9223 

~RS R8 

16.3349 


a See reference (12). 

R is assumed to be in a-units (a = 2.556 A) , ^ V(R) is obtained 
in °K . 

c k = Boltzmann's constant = 1.38062 ergs/°K . 
x = R/1.8302 , c = 6.1277 . 

e This expression is valid for large R values, the use of R = 3.1 
a-units was chosen in our calculation by convenience. 
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Table I V„2.» Radial distribution functions for 3 He and 4 He at 
T = 4°K , a 


R(a-units) 

g( 3 He) 

g( 4 He) 

R(a-units) 

g( 3 He) 

g( 4 He) 

0.700 

0.000 

0.000 

2.000 

1.251 

1.282 

0.800 

0.008 

0.006 

2.100 

1.205 

1.217 

0.900 

0.076 

0.076 

2.200 

1.164 

1.164 

1 .000 

0.284 

0.330 

2.300 

1.129 

1.123 

1.100 

0.606 

0.755 

2.400 

1.101 

1.092 

1.200 

0.931 

1.181 

2.500 

1.078 

1.069 

1.300 

1.177 

1.480 

2.600 

1.060 

1.052 

1.400 

1.326 

1.626 

2.700 

1.046 

1.040 

1.500 

1.395 

1.655 

2.800 

1.036 

1.031 

1.600 

1.408 

1.613 

2.900 

1.028 

1.024 

1.700 

1.388 

1.536 

3.000 

1.022 

1.019 

1.800 

1.348 

1.447 

3.100 

1.017 

1.015 

1.900 

1.300 

1 .360 





3 la-unit = 2.556 A . 


Figure IV. 1. Pair Distribution Functions for Helium at 4®K and 
322°K as a Function of R . 



CHAPTER V 


RESULTS AND DISCUSSION 


A. Results 

A problem worthy of attack 
Proves its worth by striking back. 

—A Grook by Piet Hein 

In Chapters II and III we have discussed the calculation of the 
polarizability functions Aa and @ for the helium diatom. The 
results of the calculation are summarized in Eqs. (3. 31 )-(3. 34) and 
in Tables III. 8 and III. 9 and Figures III.l and. III. 2. In Chapter IV 
have discussed the helium radial distribution function g(R) , which 
is given in Table IV. 2 and Figure IV. 1. We are now ready to put 
these results together to calculate the experimental quantities 
introduced in Chapter I. 

The macroscopic quantities are related to the corresponding 
molecular properties through integrals weighted with g(R) . The 
second dielectric virial coefficient (see Eqs. (1.34) and (1.39)) 
is an average of the mean incremental polarizability Aa . The 
second Kerr virial coefficient (see Eq. (1.45)) is an average 
of the anisotropy g* . The Raman depolarization ratio 6 (see Eq. 
(1.58)) reflects both Aa 2 and g 2 , but is dominated by g . Thus 
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the experimental quantities are complementary with respect to the 

components of the polarizability tensor. 

The integrals involved in B , B„ and D were evaluated 

numerically by Simpson's Rule quadrature. The final integration 

grids were chosen by varying the number cf points in the grid until 

further reduction of mesh size left the value of the integrals 

unchanged, to the number of pieces to which the results ere reported. 

The upper limit at which the integration was truncated was determined 

in the same manner. It was found sufficient to use integration steps 

of .04 a-units over the range .65 to 4.55 o-units. The B results 

e 

were checked by transforming the correspond ing range of integration 
to the interval [-1,1] , and numerical integrating the result using 
a 64-point Gausian quadrature. The conversion factors from atomic 
units to the experimental units is given in Appendix E. 

The results of our calculations cf B , , ard 6 are 

€* N 

presented in Table V.l. They are summarized as follows. Cur 

calculated value of is within the rather large experimental 

error bars at 322K. We are also in order of magnitude agreement 

with the results at 4K, but we predict the isotopic dependence to be 

opposite to that observed experimental!; ; we shall discuss this 

discrepancy in further detail below, where the quantities B U and 
L G 

will be explained. 

Our calculated at 300K verifies the null experimental 
result. We predict that the contribution of the second v’rlal term 



Table V.l. Dielectric and light scattering properties of fluid helium. 


Property 

T = 322°K 

T = 

^He " 

4°K 

Second dielectric virial 

-.094 cm 6 /mole 2 

-.025 cm 6 /mole 2 

-.030 cm 6 /mole 2 

coefficient, B e 

(-0,93 cm 6 /mo)e 2 ) a 
(-.06 ± .04 cm 6 /mole 2 )k 

(-.030 em 6 /mole 2 ) c 

(-.023 cm 6 /mole 2 ) c 


-.086 cm 6 /mole 2 

-.012 cm 6 /rnol e 2 

-.015 cm 6 /mole 2 

B e 

-.135 cm 6 /mole 2 

-.034 cm 6 /mole 2 

-.041 cm 6 /mole 2 

Second Kerr virial 
coefficient, 

2.4 x 10 -15 esu 
(3.1 x 10-15 esu) a 

1.2 x 10" 13 esu 

1,4 X 10-13 eSU 

Depolarization ratio for 
Raman scattering, D 

y 3 + -37 
V 3 + . 44 a 

ti + .08 

% + .07 
(1.33 ± . 20) d 
(1.39 ± .20) d 


Calculation of reference (27); temperature used was 300°K. 


Experimental results of reference (la), 
c Experimental results of reference (lb), 
d Experimental results of reference (36), 
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to the observed Kerr constant of helium to be ,1% at room temperature 
and pressures of 50 atm., which is not observable. However at 4K, 
the B^ is larger than the room temperature result by a factor of 
50, and the second virial term contributes 14%, at a pressure of 
1 atm. Thus it may well be possible to detect the at 4K. 

Our calculated value of D differs from the purely geometrical 
value of ^ by only 7%. Thus, it is not likely that Raman scattering 
from Aa can be observed by measuring depolarization ratios. 

We turn now to a consideration of errors in our calculation. 

In particular we wish to assess the significance of the discrepancy 
with experiment in the isotopic dependence of B e at low temperatures. 
The values of B e reported in Table V.l were obtained using the 
Hartree-Fock polarizability curve given in Eqs , (3.31) and (3.32). 

Thus, the use of our computed Act curve can be questioned on two 
main points: there is a basis set error associated with the fact 

the Hartree-Fock-Roothaan prescription is carried out with a 
(necessarily) finite basis set; and, there is an inherent (correlation) 
error associated with even the exact Hartree-Fock result. In Section B, 
the former point is discussed, and our calculation is compared to that 
of O'Brien et al_. In Section C, the correlation error is discussed, 
while Section D contains our conclusions. 

With respect to errors in g(R) , we have pointed out that at 
2°K, our g(R) values are in agreement with the calculation of 
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reference ( 7 ?) to at least three decimal places. VeHrtJor o*' 
computational parameters in the direction which would c’r.dvce a more 
accurate g(R) produced changes in, at worst, the fiftr. decfma' 
place. Further, the method used here was shewn to give good values 
of the second pressure virial coefficient in the calculations of 
reference (7A). 

B. Basis Set Error 

As has been pointed out in earlier chapters, one of the main 
objectives of our work was to compute the Hartree-Fock limit fo- the 
polarizability of the helium dimer for those litter nuclear sepa^st^ors 
which contribute to certain experimental observables. 

By definition, the Hartree-Fock limit for a giver property is 
that which is obtained through the exact solution of the Hartree-Fock 
equations. Except for field-free atoms, however, the Kertree-Fock 
equations are generally solved as we have done, via basts set 
expansion of the molecular orbitals, using the Hartree-Fock-RoGthaan 
scheme. This means that the Hartree-Fock limit corresponds to zero 
basis set error in this analytical method of solution. 

In terms of basis set selection procedures, the problem of 
obtaining the Hartree-Fock limit for field-free diatomic molecules 
is solved, and, in fact, actual calculations have been performed for 
many diatomic molecules, including the one of interest here, He 2 . 6 ’ 
From the viewpoint of perturbation theory, one should require that 
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the Hartree-Fock limit of the perturbed problem should relax to the 
field-free Hartree-Fock limit as the field is turned off. Thus, if 
the effect of the field is to couple the zero-order basis 
functions to functions of a different symmetry (this is the case for 
fields applied perpendicular to the internuclear axis of a diatomic 
molecule), it is clear that this can be accomplished by augmenting 
the optimized field-free basis set with functions which describe 
their distortion by the applied field. For cases in which the 
polarization functions are of the same symmetry (as is the case for 
fields applied parallel to the internuclear axis of a diatomic 
molecule) one follows the same prescription. The reason for this is 
that if the field-free wavefunction truly represents the Hartree- 
Fock limit, the effect of the polarization functions on the energy 
and properties for the zero-field case will be minimal. 

The manner in which we selected the basis set in this work, was 
described earlier (Section C of Chapter III). With respect to 
assessing the success of our approach in light of the foregoing 
discussion, the following three points are significant: 

1. Our zero-order set was essentially a Hartree-Fock limit 
set, with orbital exponents optimized (by Kestner) at each 
internuclear separation. 

2. Any reduction in either the zero-order or polarization set 
was made only after it was justified by numerical 
experimentation. 
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3. The Hartree-Fock limit of the polarizability was obtained 
for both the separated and united atom limits, using atomic 
extensions of the diatomic prescription. 

Based on his own basis set variations, and on the results of 
the more elaborate calculations of Gilbert and Wahl, Kestner asserts 
that the total energies computed with his wavefunction differ from 
the Hartree-Fock limit values by less than one part in 10 7 . To this 
extent, one can be confident that the choice of Kestner's wavefunction 
is a reasonable one. Although we used the same basis functions for 
both the tjg and occupied orbitals, whereas Kestner used 
different ones for each, the effect of this difference is expected 
to be small, since the corresponding change in the field-free energy 
is on the order of a few parts in 1G 6 . This means that as the field 
is turned off, we still relax to a function which obtains essentially 
the Hartree-Fock limit for the energy. 

It, of course, does not follow that any function which yields 
the Hartree-Fock limit for the zero-order energy is as good for 
second-order property calculations as any other which also yields 
this limit. However, additional credence is lent to our results 
because our basis goes smoothly into that for the separated atom 
limit, for which we obtain the correct coupled Hartree-Fock limit 
polarizabil ity. 

In comparing our results to those of O'Brien et al_., three points 
are especially noteworthy: first, in their calculation, they used 
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the same (Gaussian) basis set for all values of the internuclear 
separation; second, the difference between the Hartree-fock limit 
for the total energy and their total energy is on the order of a 
part in 10 1 *; third, they failed to obtain the coupled Hartree-Fock 
limit for the polarizability of the helium atom. 

The use of a single basis set for all internuclear separations 
has the advantage that smoothness is automatically built into the 
polarizability curve. For the case in which the single basis is an 
atomic set, one thereby introduces the so-called distortion error. 61 
Using our basis set, this amounts to an error of about .003 in 
ct^ at R - 6.0 a 0 . This and the use of a Gaussian basis set are 
probably the major source of error in the computed, zero-order 
potential curve of O'Brien et al . 

The most significant difference between our results and theirs 
is that they do not obtain the Hartree-Fock limit for the polarizability 
at the separated atom limits. In fact, while our values of a (R) are 
slightly below theirs, our values of Act(R) are about the same; this 
is because their value of 2a 0 is greater than the Hartree-Fock 
limit by about the same amount that their values of a(R) are 
greater than ours. 

We estimate that our results are within about .002 of the 
Hartree-Fock limit of the polarizability for R values between 4.0 a 0 
and 5.5 a„ , and are within about .001 a 0 for R values between 
5.5 a Q and 6.0 a Q . Curves for Aa which correspond to our error 
estimates are given as follows: 1 
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a U (H) « -14.5expM.48R) + .C085/R , R < 6.0 a s (5.1) 

a L (R) ^ -12.9expM.45R) - ,0127exp(- ,462R) , R < 7.0 a .(5.2) 

For large values of R we use the po^nt dipole result, Eq, (1.37) 
and (1,38) , These curves are shown as the bounding curves in Figure- 
Ill. 1. The corresponding values of the second dielectric virial 
coefficient are denoted by B U and B L in Table V.l. In this 

E C 

table the "U" superscript refers to the result obtained using Eq. 
(5.1) while the "L" superscript refers to the result obtained using 
Eq. (5.2). We should emphasize that these- bounds are based on "best 
guesses" of errors, and are not quantitative. The anisotropy curve 
which is of interest for the calculations is giver in Figure 
III. 2. 

It is apparent from Table V.l that the magnitude of B £ depends 
markedly on which Aa(R) curve one uses; in fact, the value of B L 
at room temperature is seen to be out of the experimental range. It 
is equally apparent, however, that the isotopic ordering of the low 
temperature results is the same regardless of which Aa(R) curve is 
used. Thus, we conclude that basis set error does not account for 
the low temperature isotopic discrepancy. 
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C. Electron Correlation Effects 

This section is not intended to be a general comment on the 
effects of electron correlation on molecular polarizabilities, but 
rather, it concerns itself with whether or not this effect can 
account for the isotopic reversal of our values of the second 
dielectric virial coefficient for 3 He and 4 He relative to experiment. 

To get a better feeling for the problem, one can examine 

Figure IV. 1, which shows a plot of the g(R) values we obtain for 

3 He and 4 He at 4°K, along with that for room temperature helium gas. 

Our Aa(R) points are plotted in Figure 111.1. The isotopic' differences 

between the g(R) curves at this temperature is due to mass effects: 

that is, the larger mass of 4 He results in an increase of 4 He density 

over the well, while the smaller mass of 3 He allows it to tunnel 

further into the barrier than 4 He. The order of the B values 

e 

which results is then a resolution of the question whether or not 
the slightly greater penetration into the barrier by 3 He to regions 
of more negative Act(R) can overcome the significant increase in 
4 He density over the well. 

In order to effect the reversal of our results by invoking 
correlation effects, one must argue that the effect of correlation is 
to "pull up" the Aa(R) curve for internuclear separations over the 
well and "pull it down" in the tunneling region. This latter effect 
is limited for two reasons: first, Aa(R) is bounded from below 

because a(R) must be positive for a molecule in its ground state; 
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second, even within the constraint introduced by the above point, a 
decrease in the Hartree-Fock Aa(R) curve in the tunneling region 
for a broad range of R values will destroy the agreement between 
the room temperature experiments and the Hartree-Fock calculation. 
Thus, in order to agree with experiment, one would have to require 
that the effect of correlation on the incremental polarizability 
curve is to make Aa(R) positive over the well, negative and 
sharply decreasing in the region where g{ 3 He) > g( 4 He) , and 
negative but increasing in the region where the classical g(R) is 
large, and the quantal g(R) curves are small. 

The results of earlier chapters indicate that the correlation 
contribution to the mean incremental polarizabil ity in the dispersion 
region is positive, and is given by: 

a Corr (long ~ range) = < 38 - 98 ' 9.24)R“ 6 aj . 

Here, 38.98 a® is the accurate A 6 result of Chapter II, and 
9.24 a Q is the corresponding point dipole value. If one assumes 
that the multipole expansion is valid into R = 6.5 a 0 , this 
results in slight increase in B £ for 4 He relative to 3 He. 

For separations in the overlap region, the correlation 
correction is not known, and any conclusions are therefore 
speculative. To obtain an idea of what one is up against in this 
region, we denote the exact (correlated) result by a , and the 
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Hartree-Fock result by a , Act (correlated) by Aa , Aa(Hartree- 
Fock) by ta, , and consider the following: 

Aa a(R) - 2ag > 

Aa - a(R). - 2ag , 


or. 


Aa = Aa - .1223^ f (a(R) - a(R)) , 

where we have used the known result, a 0 - a 0 - 0,061a G . GS The 
sign of Aa in a given range of R is determined by the sign of 
the last term on the right-hand side of this equation, since the 
first two terms are both negative. Intuitively, one might expect 
this term to be negative in the region of the van der Waals well, 
corresponding to the greater diffuseness of the Hartree-Fock 
wavefunction compared to the correlated function. This is the case 
for most atoms and for H 2 at its equilibrium separation. However, 
definitive statements regarding He 2 can be made only by actually 
performing the correlated polarizability calculation. 

In an attempt to match the isotopic dependence and magnitude of 
the experimental data at 4°K and 322°K, we adopted several ad hoc 
models. These models served to distort the Hartree-Fock Aa(R) 
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curve in a manner which is tern •• stent with the ideas o'.stussed 
earlier in this sect ten.. T he it irg _.»(R) cu've »equ1red a 
maximum of JO al in the -ag-.cn et the «e^, and a minimum of about 
-„30 a^ in the tunneling 'eg on while this general shape is not 
implausible, the actual numbers n velvet, •■equ'^e that, tne correlation 
correction to the mean ™c'ement<J pc ? a'- ?ab ?T ' : ty be an orde v of 
magnitude larger than a.>f Hardee- Feck} 

The effect of correlation of the mean incremental polarizability 
on the properties other than B lasted r n Table VJ, 1$ minimal, 
so that no insight as to > ts magnitude can be inferred from those 
experiments. 


D . Summary and Conclusions 

The largest parr of the d’M.v* v-cn has dealt with the dielectric 
constant experiments, as thts is the only plate where a conflict 
between the present work and experiment, has appealed To summarize 
our dielectric constant results, we see that we get good agreement 
between our computed results and those of the room temperature 
measurements, and order of magnitude agreement with the low temperature 
experiments, but with reversed < so topic dependence. Ad hoc 
modification of the incremental mean peParizabil it.y curve within what 
we take to be reasonable limits dees not reverse this result. Our 
estimates of the other properties which we considered all agree with 
experiment where such experiments have been performed; agreement with 
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other calculations of the high temperature properties is also good, 
although the corresponding polarizability curves are somewhat 
different* In this respect, we believe our calculations represent 
the Hartree-Fock limit more closely than the others do* 

Further, it is our opinion that the disagreement with experiment 

is a consequence of using too small error assignments in the difficult 

experiments. This reflects on the analysis of the experiment in two 

ways: first, in the region of low densities where the virial 

expansion is certainly valid, the data exhibits a large amount of 

scatter, so as to render these points relatively useless as data in 

a linear least-squares fit* On the other hand, in the region in 

which the scatter is smallest the number densities are quite high. 

At these densities, a two-term virial expansion of the equation-of- 

state is known to be a poor approximation* Thus the density 

coefficient quoted by Kerr and Sherman may be an average slope, as 

distinct from the initial slope which defines B * Neither of 

£ 

these problems appear in the room temperature data of Orcutt and 
Cole. 

Our results indicate that a priori calculations of the molecular 
properties of helium can now make meaningful contributions to the 
interpretation of experimental data, and we believe that renewed 
experimental and computational efforts on the dielectric properties 
of helium is very desireable* In this connection, additional 
calculations on X-ray structure factors for liquid helium; third 
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pressure viri'al coefficients for the gas, and Raman scattering line 
shapes are presently being carried out in this laboratory. 



APPENDIX A 


SYMMETRY REDUCTION PROCEDURE USED 
IN THE LONG-RANGE CALCULATION 

The reduction of Eq. (2.32) to (2.48) is long but 
straightforward. We shall illustrate the technique for a simple 
example, and then present a table of results. Consider the single- 
center matrix element. 


I = <4>oly m v_ m Uo > » (A.l ) 

where y ±m ~ Y^ ±in , and Y £m represents the usual spherical 
harmonic. In this work, we used the phase convention, 
y* = ) m y 

The intent of the symmetry reduction is to relate the matrix 
elements of the form (A.l) which have m t 0 to those which have 
m = 0 , so that Eq. (2.32) can be simplified. To effect this, it is 
natural to make use of the Wigner-Eckart theorem. 45 This theorem 
states that, in a standard representation {J 2 ,J } , whose basis 
functions are denoted by |t,J,M> , the matrix element 
< tJM|T^ K) |t'J-M'> of the q-th standard component of the k-th 
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( Kl 

order Irreducible tensor operator, T v ' , is equal to the product 
of a Clebsch-Gordon coefficient, < J ' KM 1 q 1 0M> and a quantity which 
is independent of M, M', and q . That is, 

<x]MlT^ K} |t l J , M , > = <Tj||T (K) |!T'J l ><J , KM'q|JM> , (A. 2) 

H 

( K) 

where the quantity <tJJ|T v 7 1 [ -r * J ‘ > is known as the reduced matrix 
element, and can be evaluated by using a particular value of M for 
which the left-hand side of Eq. (A. 2) can be easily evaluated directly. 

In the following, we ignore the radial part of the functions y m , 
as they are independent of m , so that we take, 

~ (A. 3) 


Then, 


* Mi“<Y lfa *o|Y lt .*0> (A.4) 


(-l) m <Y 1 


,m y O' Y 0,ol Y i, m 


(A. 5) 


and 


<( * , ol Y i J o Yi,oU 0 > = ^ <Y i,o ♦ol Y o,ol Y i,o 4>o > • (A. 6) 
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Now, since the 2k+l spherical harmonics (q - -k, -k+1, 

k-1, k) considered as operators, are the standard components of an 
irreducible tensor operator of order k , , we can write 

Eq. (A. 5) as: 


(-l) m ^7r <Yj 0 J Y 0 , 0 I Y ! m < Po > - (-1)“ rfir <Y^ 


m 


T °»°l Y l,m 4, 0 > 


(A.7) 


Further, since <p 0 has S-type symmetry, the right-hand side of this 
last equation is of the form of Eq. (A. 2), so that we can apply the 
Wigner-Eckart formula to Eq. (A.7) for the two cases m =•■ o,l : 


<y i,i4>oI t o,o( y i 1 i4 > o > 


' 1 0 l' 
~1 0 1 

/ 


<T »1 1 |T 0 I I T ' »1> , 


(A. 8) 


and 


<Y i,o | l ) ol T o,o|Yi,o < l > o > 


(- 1 ) 


1 0 V 

0 0 Oj 


< T » 1 i f To I l T ' *1> 


(A. 9) 


where we have replaced the Clebsch-Gordon coefficient of Eq. (A. 2) 
with the corresponding Wigner 3-j symbol. Thus, 




( - 1 ) /4ir 


1 ° 11 

-101 <Xl1i|T o | l T '* 1> * 

k- 4 i 
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and 


1 0 I, 


♦ olYi.oYi.oko> " H)/?« <t,1| | T 0 1 ( t’ . 1> . (A, 11) 

1“' 0 lj 


but since 


and 


we get: 


fl 0 1] 
0 0 o! 


/3 


i 1 ° lj , 
:-1 0 1, /3 


<^0 I 


> - 


Iwiv-i l4>o > 


(A. 12) 


We needed a number of relations such as Eq. (A .12) . They are 
collected together In Table A, These results can also be obtained 
by using explicit relations for the 's or by raising-lowering 
operator techniques. 

Let us now outline the symmetry reduction of Eq. (2.32) for 
A 6 . First we note two simplifications: 

1. The dipole moment operator can be replaced by y a , the 
dipole operator for atom a , since the cross terms involving p a 
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and y b vanish* The final expression is symmetrized by multiplying 
by (1 + P , where P ^ exchanges the labels of a , b » 

2. The resolvent r 0 is spherically symmetric about each 
center, and hence the excited atomic states behave like the 
spherical harmonics Y Am . 

The dipole-dipole interaction operator can be written 

i; _ n.ab 1 , a b 0 a b 

*3 " %V + P_ " P + - 2p 0 y Q , 

where y 0 ~ v z » = -u x - iy y and y_ = y x - iu y . Then one of 

the nonzero terms in A 6 is 

<< f > 0 I ^3 r 0 V 3 r 0 P+ r 0 y_Uo > = O+P ab ) < ^0 [ V 3 rS b V 3 Vq y® rg y*|* 0 > 

- KO +p a b) <( *oli J + v b _ r o b K P+ r o p+ rS uf f + • ... 

For the excitations involving atom b , we have 

^ l-iYiixfYn) Y n > - -<Y 1 . 1 Y 11 > 2 = -<Y 10 > 2 . 

For atom a , the middle resolvent can involve either excited 3(Y 00 ) 
°r D(Y 2m ) states. Thus, there are two possibilities: 
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<Y ii Y i-i >< ^i-i Y i-i Y oo ><Y oo Y ii Y i"i >< ^i-iYi-i> 


= (-1 ) (+1 ) (-1 ) {+1 )<Y iq y io ><y xo y io y oo ><y oo y io y i o >< Y io Y io 3 ” 


and 


<Y U Y X - 1 ><Y 1- 1 Y 1- 1 Y 20 ><Y 20 Y 11 Y 1-1 ><Y 1-1 Y 1-.i. > 

("1 ) (-/£) ( ) <Y 10 Y 10 ><Y 10 Y 10 Y 20 ><Y 20 Y 10 Y 10 ><Y 1 0 Y 1 0 ^ • 

This gives 



^olvo Vo r o y o vg r o 4 rS woUo > s 

+ H- << t > 0 1 Vo Vo r 0 Vo Vo Vq pg Tq Pq|(J)q>p , 

where the subscript S and D refer to the symmetry of the excited 
states in the middle resolvent. 

In the way we have just illustrated, we reduced Eq. (2.32) to 
Eq. (2.48), with the coefficients given in Table II. 1. 
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Table A.l. Symmetry reduction of sperical harmonics. The angular 
brackets < > denote spherical averages. 

< YhYi_i> = -<Yi 0 > 

< YioYuY 2 -i> - - ~2 

< YnYi_iY2o > = 1 2 < YioY20 > 

= -<y| 0 > 


<Y 2 1 Y 2- 1 > 


APPENDIX B 


DERIVATION OF INTEGRALS REQUIRED IN THE 
LONG-RANGE HELIUM CALCULATIONS 


The basic reference for integral evaluation for the long-range 
calculations is the atomic integral work of Calais and Lowdin. The 
general integral for which they derive closed-form expressions is: 


Uc|bd) = ff f a ( r i>\ imi ( f >i’U)f b ('- 2 )y^ ni2 (e 2 ,^)h(r 12 ) 
x f c (r ‘ )y i,m, (6 i-»i) f d< r 2) Y umu (6 2-*2 )d !:i d t2 


(B.l) 


In reference (51) this integral is reduced to a sum of simpler 
integrals of the form: 


r mti 

U 


ff 


f(ri)Y 


£m 


(e 1 ,4>x)h(r 12 )g(r 2 )Y Ay {e 2 ,* 2 )dr 1 dr 2 . (B.2) 


For our purposes, the functions f,g, and h are of the form r v e' yr 
We use the following notation. 



(B.3) 
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where 

[0,M] *• (a+b)(a+c)(b+c) • < B - 4 ) 

and 

aiK 35T = [k,£,m] . (B.5) 

Thus, knowing the fundamental integral [0,0,0] , one can generate 
formulas for the set [k,t,m] by carrying out the differentiations 
indicated in Eq. (B.5), and (ac|bd) can be computed. The general 
expression for the [k,A,m] is given in reference (51). 

In the present calculation, construction of the Hamiltonian 
matrix requires, in addition to the , expressions for the basic 
integral : 

~ | ^ r i^i,m^i ’^i) e ^ 9( 1 "2)Y^(02» ( f’2)h(r l 2)dri dr 2 . 

(B.6) 

One purpose of this appendix is to derive a computationally useful 
expression for this integral. 

Following Calais and Lowdin, we first integrate over $ 2 , 
and ej and express the remainder as an integral over r^, r 2 , and 
the relative angle, e 12 = e 2 - e x . To accomplish this, we first 
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rotate the coordinate system so that the polar axis for r 2 coincides 
with rj . Denoting the new polar angles as e 12 and $ 12 , we have 
that Y X(j (e 2 ,<}> 2 ) becomes: 


Y Ay ( 0 2 »<{ | 2) 


■ l e 




m 






(B.7) 


Here, the factor e 1 ^ 1 d* , ( e i ) is an element of the rotation 
matrix representing rotation through the Euler angles a = » 

3 — 6 i » y = 0 . Using Eq. (B.7) in (B.6), we can integrate over 
♦ 12 and see that only the m' » 0 term contributes to the above 
sum. Thus, since, 

V e -*> ■ Vcoseje 1 "* . (B.8) 


and 



we get 


(B.9) 


/ , sm+u p(a-u)! 2A+1 (2t+l ) U-m) ! -, l /o 

L TX+5yi 4^ ^U+m)! - 1 

x Jf(ri)g{r 2 )h(r 12 )P £m (cose 1 )P Xm (cose 1 )P x (cose 12 )e i(m+n+lj) ^i dt , 

(B.10) 


mpn 

y U 
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where 


dt * r* r* sine 1 sine 12 dr : dr 2 dr J2 de l de 12 d^ . 


Integrating over yields a factor 6 ^ 9 so that y^t can 

m-ry , -n J6A 

be written, 


myn 

hx 


“ 2 " S m+p,-n f f < r ll r l dr l f S< r 2> r 2 ir 

* ^ J*l 


f ir 


h(ri2)P >i (cose 1 2)sin0 1 2 de 12 


(B.ll) 


f- 1 )" <2 4t(Um)T ?! J 2 " f P *,n<‘'>V u > «" • 


Comparing this result to Eq. (17) of reference (51), we get 


myn _ 
Y ZX 


• 2 ’ (‘ ■ , J“>V'> d “] 


x r 


m+n,y 

ZX 


(B. 12) 


Hence, y^ n can be done by computing v via the method of 

Calais and Lowdin, and modifying it by the above multiplicative 
factor. The integrations involving the associated Legendre functions 
were done explicitly for the particular P 's which we required. 
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This result, Eq 0 (B.12), can be duplicated by expanding the 
factor in terms of spherical harmonics, and using the 

addition theorem for spherical harmonics to reduce the result to the 
form of r^ n,y . 

As mentioned earlier in this appendix, the integrals of Eq. (B.6) 
can be reduced to sums of integrals of the form of Eq, (B, 3), with 
Eqs, (B.4) and (B,5) forming the basic computational units. For our 
purposes, however, we also need the additional integral: 


I 


r I 1 e 


-kr 


a 

1 dr x 

r 


e- kr 2 dr 2 




r l +r 2 

I r l -r 2 I 


r i2 dr 


12 > 
(B.13) 


for -1 £ b ^ s and 1 £ c £ t . 

We proceed by integrating over r X2 first, then over r 2 , and 
finally over r x , using standard expressions for the intermediate 
integrals which occur. The final result which is obtained for the 
integral I is: 


= Ii + I2 » 


I 
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where 


- { ( c4i)k^K ci4( ' i)t4?i(bH4,)i *" 2 + i ci+( - i)p+1] 

p - x 

, IctULlbielUc^llt, . 1 ] .. f Y 1 


p!(c+l-pTT 


p=l q -- 0 


[c»1)i (bt p)l (bH -q)j 1 
pT(c+l-pyi~Tb+p-q) I 2b+c.-q+rj ' 


(B.14) 


For I 2 , we have different results depending on whether 
b - -1 , or b -] i 


bj 

h 


-1 


" H+T FW l f {-1) C+2 ]*n2 + kSTF+2 * i [!+(-'< ) P+1 ] 

'P“l 

* pifiFr^irr o>* c+ i- P ):n - - T F r 

P“1 q=0 


(B.15) 
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b > -1 : 

>2 * ciTpIcTZ'fC 1 ' jjfjTjMclCH-l) 1 *" 2 ] + f [l+(-') P+1 ] 

( r xi \ 8 i c+1 p-2 

x pTt §^{-p)l (P-1)» (b+c+l-p) ! [1 - ^rp] -IV (B.16) 

p=l q=0 


X 


[> + <- 1 > p+ 1 V& 


(p-1 ) i ( b+c-q) ll 
(p-q-lj! 2 b+c+ l-qj ‘ 


Eqs. (BJ2), (B.13), (B.14), (B„I5), and (B..16) taken together 
with the results of Calais and Lowdin comprise all the necessary 
expressions for computing the integrals needed in the long-range 
helium calculation. 



APPENDIX C 


THE BISON COMPUTING SYSTEM 


BISON is a self-contained FORTRAN computer system for the 
calculation of analytic self-consistent-field wavef unctions, 
properties, and charge densities for diatomic molecules, using a 
basis of Slater type orbitals* BISON was written by A* C* Wahl, 

P* J, Bertoncini, K„ Kaiser, and R„ H* Land of Argonne National 
Laboratory and is described in reference (64)* In the initial 
construction of this package, care was taken to guarantee that the 
code would be computationally optimized, modular, machine 
transferable, and capable of continuous growth and revision* For 
these reasons it is particularly attractive for our purposes* 

Before the program could be modified to allow for computation 
of diatomic polarizabilities, the standard working version had to be 
made operative on the University of Wisconsin UNIVAC 1108* The 
version of BISON with which we began, consisted of ten overlays 
(independent program segments) which are defined as follows: 

roain - basic executor, remains core-resident throughout 

the computations; 

overlay 1 - input driver and data-screening segment; 
overlay 2 - compute J-type exchange integrals; 
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overlay 3 - compute J-type coulomb-hybrid integrals; 
overlay 4 - compute K-type exchange integrals; 
overlay 5 - compute K-type coulomb -hybrid integrals; 
overlay 6 - compute one-electron integrals; 
overlay 7 - construct P and Q matrices; 
overlay 8 - carry out SCF iterations; 
overlay 9 - compute properties; 
overlay 10 - contour hunting and drawing. 

A description of the notation and definitions of terms can be found 
in the BISON manual. 

The initial problem was simply enough storage space to store 
intermediate results. The data bank (d-bank) size of the above 
overlays was on the order 130 thousand machine words (130k) 10 (double 
precision version) except for overlays 1, 7, 9, and 10 which were 
considerably smaller. The data divided up among several arrays, 
with one array of dominant size in the two-electron integral overlays 
and several arrays of similar size in the one-electron integral 
overlay. In each case, however, the sizes were such that no matter 
how the arrays were arranged in core, at least one array would begin 
beyond the nominal memory size (65k) ig » This comprised a main 
problem in adapting the program to the 1108; because the UNI VAC 1110 
series of computers uses a 16-bit address field in its instruction 
words, thus the base address of any array must be less than (65k) 10 
relative to the beginning of the program. However, special data 
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arrangement, procedures can be used to extend this limit to (256k ) kj = 
This is done by collecting a sufficient number of arrays into a 
master array using the FORTRAN V DEFINE procedure, such that the base 
address of this master array is less that (65k)io . 79 The words 
beyond the (65k) 10 boundary are addressed by the operating system 
via the base address plus an index offset. This clearly implies 
that this master array must be stored as the (physically) last part 
of the d-bank for the given overlay. This is accomplished by 
declaring the master array to be in labelled COMMON, and manipulating 
the position of this COMMON block in core using the EXEC-8 MAP 
processor (MAP is the local name for the EXEC-8 COLLECTOR). 80 

These points are best illustrated by an example from the 1108 
version of BISON. Figure (C.l) is a partial listing of the source 
code for subroutine SUT, which is the main subroutine in the 
one-electron integral overlay. Figure (C.2) is a listing of the 
output of the MAP processor for the same overlay. 

In Figure (C.l), the COMMON block of interest is EXCORE (a); 
since this is a double precision program, the array STRAGE requires 
(130,000) 10 words. Part (b) of Figure (C.l) shows the actual arrays 
which are contained in STRAGE, and the manner in which they are 
imbedded. The statement which defines the function MM(I,J) is of 
additional interest, because it illustrates the procedure used to 
imbed a two-dimensional array in a one-dimensional array. The 
offsets given in Part (c), are not to be confused with the operating 
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Table C.l. Illustration of the Use of the DEFINE Procedure in 
BISON. 
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SEGMENT 5UT * Q 1 4070 0177*10 05R761 502377 



FOLLOWS “SEGMENT 

MAIN 



« 

OFFSET 

(COMMON 

BLOCK) 




054761 

05R763 

SUTMAT 

(COMMON 

BLOCK) 




05R76R 

1Q0RQ3 

SUTPUT 

(COMMON 

BLOCK) 




100RQR 

100SR7 

SUTDUM 

(COMMON 

BLOCK) 




100550 

100713 

“WfDAT 

(COMMON 

BLOCK) 




1007 14 

102077 

SUTCOM 

(COMMON 

BLOCK) 




102100 

103537 

SUT 



1 

014070 015662 

0 

1035R0 

103770 




3 

INPUT 

2 

FLOAT 





5 

FLAGS 

R 

PARM 





7 

SUTDUM 

6 

SUTOAT 





? 

EXCORE 

8 

SUTPUT 





11 

SUTCOM 

10 

01PLM 





13 

OFFSET 

12 

suthat 







IN 

1 BUCKS 

; • 

CXTTSS - 



1 

1 015663 016175 

0 

103771 

104113 




3 

SUTOAT 

2 

FLOAT 





““ 5” 

— SUTCOM 

- B. M 


■■■ ' 

STOSUT 



1 

016176 016535 

■El 






3 

SUTPUT 

2 


m 




6 

SUTCOM 

4 

EXCORE 







6 

OFFSET 


LEGSUT 



1 

016536 017354 

0 

10R176 

10R323 




3 

“EXCORE 

2 

FLOAT 

; 






R 

offset 


ADDSTQ 



1 

017355 017111 

0 

10R32R 

104335 






2 

SUTDUM 

3 

blkctc 





0 

10R336 

104336 






2 

mmnjM 


DPARAM 



1 

017412 017660 

P 

10R337 

10RR3S 




3 

D1PLM 

* 

EXCORE 





5 

OFFSET 

R 

SUTCOM 


PUTLEG 



1 

017661 017740 

0 

10RR36 

10RR57 




3 

SUTPUT 

2 

PARM 

* 

FxtORE 

(COMMON 

“BLOCK)' 




10RR40 

502377 



\ 

_\ 


-t- 
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Table C.2. Illustration of the Use of the MAP Processor to 
Position Data. 
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SUBROUTINE SUT 

IMPLICIT REAl*8 (A-H.b-Z) 

RE A L* 8 MM 

R £ A L • 8 M A TEL : 

REAL GRID 

COMMON /FLO AT /FLINT! 1 7 6 ) .FACTL ( 57 j 

COMMON/ INPUT/ Z 11 00 ) . NN ( 1 00 ) , NL U 00 > , NM ( i 00 ) . NNEU # NS I MP , NX I .NETA^ 

1 NSYM.NBIAS, ICH.NLAMDAUO) 

COMMON/PARM/ Z A , ZB , RR , RM A X , B ET A , MMM I N , MMM AX , LLM I N , LLM AX , I PM AX , NMAT 

~c omo^/ra g s^^ s tq 7 tfun~ctt 

COMMON/SUTDAT/GRIDU20) 

common/sutduH/ toumhoo) 

COMMON/SUTPUT/ jP LACE! 100) 

COMMON/EXCORE/STR AGE! ASOQO) V g. ' " 

COMMON/OIPLM/ I DIPL. FIELDS 

CO M MO N / S U T COM/GDPT (2.T66) , GDWT ( 2, 100) 

COMMON/SUTMAT/ M ATEL ( 1 000 . S ) _ _ 

COMMON/OFFSET/ 10FST 1 , I 0F5T2 . I OF ST 3 

COMMON/ 1BLCKS/NBLCKS 

DEFINE PLMNX ( I ) cSTRAGE ( I ) 

DEFINE CHI l n«STRAGE( I + I0FST1 ) 

DEFINE RT I)b S T RAGET I + I 0FST2 ) “ 

DEFINE COS! 1 1 • ST RAGE l I ♦ I OFST 3) 

DEFINE MM! I . J ) “ STR AGE ( J* NROWS-NROWS + l + I OFST S ) 

DEFINE ZVALU! I . J)«STRAGE< J*NRSP-NRSP+ I + 1 0FST5) 

DEFINE X V ALU I I) ■S TRAGE ( I ♦ I OFST A ) 

IOFST1«10000 ' 7 ~ 7T \ 

I OF ST 2*20000 " ~ * 

I0FST3«2500D 

1 0FSTH-30000 

IOFSTSsSSOOO | C 

I OFST & ■60000 

NROWS-7 

NRSP-1000 - - ‘ 
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system offset indices mentioned above, which are transparent to the 
user. The program offsets (c) are user defined and are equal to the 
total number of words in the individual "sub-arrays". 

In Figure (C.2), we see that the COMMON BLOCK EXCORE has been 
forced to the end of the d-bank of the overlay, SUT; this, in fact, 
forces it to the absolute end of the overlay, as the I-bank is stored 
first. This positioning is accomplished by placing the card, IN 
EXCORE (left-justified in column 1) as the last card in the set 
which defines SEGMENT SUT. Note that since EXCORE begins at (64772) 8 
and ends at ( 502377) 8 (relative addresses), it could not have fit 
into the addressable memory in any other way. 

Besides the core requirements, BISON uses eleven periferal 
files. These may be tape files, disk files, drum files, or any 
combination thereof. For production runs, the program was copied 
from tape on which the relocatable elements were stored to a 
temporary drum file, and collected. The rest of the periferals used, 
consisted of an integral save tape, and eight (scratch) drum files. 

The logical unit numbers assigned to the scratch files were: 

40, 41, 42, 49, 51, 52, 53, 54 . 

Of these, all but unit 52 were allotted a maximum of 128 tracks 
(1 track = ( 1 ,792) 10 words); unit 52 required somewhat more, the 
actual number depended on the basis set used. LUN 43 was assigned 



to the integral save tape; the 7-channel tape was taken to be of 
medium density (556 bpi). We used 7-channel tapes in order to be 
compatible with the 11 08 1 s 36-bit word structure. 

Once the original version of BISON was working, the primary 
modifications which were necessary to enable us to compute 
polarizabilities took place mainly in the integral overlays. The 
original BISON was designed to take advantage of the or symmetry 
of diatomic molecules. The symmetry is lowered when an externa] field 
is present in the perpendicular direction, and this required 
modifying some of the algorithms of integral evaluation. These 
changes consisted of four Sasic steps: 

1. Compute the dipole moment matrix elements in the one- 
electron segment. This is equivalent to constructing p*F 
for unit field strength. The basic code was extracted from 
the properties segment. 

2. Modify the two-electron integral overlays to allow for 
Slater orbitals of different m-values (including those 
with m-values of opposite sign) to be handled in a single 
symmetry block. This included insertion of new selection 
rules and phase factors. It was also necessary to restrict 
the MCYCLE 00-loop parameters of subroutines K.EX and KCH 

to a single value, as it was no longer necessary to compute 
both the K - and K -type integrals separately; 
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3. Modify the PQ overlay to take the loss of permutational 
symmetry in the two-electron integrals into account; 

4. Modify the SCF part of the program to include the dipole 
moment terms in the one-electron contributions to the Fock 
operator, and to allow for field variation in those terms. 

Various and sundry additional changes in other parts of the code were 
also necessary, but these will be discussed elsewhere. A description 
of the (free-format) input to the polarizability load module is given 
in Appendix D. 

The final verison of the program which we used, typically 
required about 41 minutes of CPU time to compute a value of or 
for a given value of R , with an 18- term basis set, using three 
different field strengths. Of the 41 minutes, 39 minutes were spent 
In evaluating integrals, with the remainder spent in the SCF 
calculation. A significant cost factor was the core usage and I/O 
changes, as illustrated in the typical breakdown shown below: 


CPU time $99.00 
10 requests 16.00 
10 word transferred 2.50 
CORE usage 54.00 


Shortly after the computations reported in this thesis were 
completed, the University of Wisconsin Academic Computing Center 
upgraded their 1108 to an 1110. Although they made the claim that 
programs which operated on the 1108 would operate on the 1110 without 
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modifications, our initial attempts to run BISON were unsuccessful 
due to decreased available core size on the 1110, However, since 
the final form of the 1110 operating system was not available, no 
attempt was made to re-organize the BISON code to fit into the 1110 
core. 



APPENDIX D 


SAMPLE INPUT FOR BISON POLARIZABILITY LOAD MODULE 


A sample of the (format- free) input to the BISON load module 
which is used in the polarizability calculations is given on the 
Figure D.l. The input, except for the second and third cards, is 
identical to that described in the BISON manual, reference (64). 

The data which is entered on the second and third cards is now 
described. 

The ninth (MASSA) and tenth (MASSB) entries on the second card 
represent the masses of centers A and B respectively, in atomic 
mass units. The eleventh (IDIPL) and fourteenth (NBLCKS) can be 
considered together as follows: 

IDIPL = 2 , NBLCKS = 1 ; computes the matrix elements of z in 

the one-electron part of the program, 
for unblocked input. 

IDIPL = 2 , NBLCKS f 1 ; computes the matrix elements of x in 

the one-electron part, for unblocked 
input. 

IDIPL a 1 , NBLCKS f 1 ; computes matrix elements of z for 

blocked input. 

160 
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The twelfth entry (FIELDS) is the field strength used in the one- 
electron part of the program, to compute the yF terms. The 
thirteenth entry (DEPTOL) is the linear independence tolerance, and 
execution is terminated if the smallest eigenvalue of the overlap 
matrix is less than this threshold. 

The first entry on the third card (fiFIELD) indicates the number 
of field strengths to be employed in the current run, and the next 
NFIELD entries on this card are the specific values of the field 
strengths which will be used. The field variation is done in the 
SCF part of the program, so that the use of the field strengths 
given on the third card presupposes that a value of FIELDS * 1.0 
is used on the second card. 
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Table D.l. Sample Input for BISON for Calculation of a ( for He 2 . 
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HE? CALCN. OF ALPHA-PERP. AT Rr6.0. BASIS OF 18 FNS. 

6.0 30 SO 50 0 0 ?. 0 ?. 0 4.0 4.0 ? 1.0 .10-6 0 

1 .001 

14 60 ?0 

0 1 07766700000 
18 ? ? 0 

2.0 2,0 

100 1.460 

200 2.641 

200 1.773 

2 1 1 1.450 

2 1 -1 1.450 

3 1 1 1.450 

3 1 -1 1.450 

3 1 1 1.655 

3 1 -1 1.655 

1 0 0 -1.450 

700 -2.641 

200 -1.723 

2 1 1 -1.450 

7 1 -1 -1.450 

3 1 1 -1.460 

3 1 -1 -1.450 

3 1 1 -1.655 

3 1 -1 -1.655 

-.0400 

*0800 
.1900 
-.0001 
.0001 
-. 000 ? 

, 000 ? 

.0001 
-.0001 
- ,6400 
.0800 
.1900 
-.0001 
.0001 
-. 000 ? 

,000? 

.0001 
-.0001 
.9700 
-.0800 
-.1900 
.0001 
-.0001 
.0003 
-.000? 

“.OOOl 
.0001 
-.9700 
.0800 
.1900 
-.0001 
.0001 
-.0003 
.0003 
*0001 
-.0001 

END OF^RUM 



APPENDIX E 


CONVERSION FACTORS FOR B AND B. 

C l\ 

In our work, we have computed the polarizabilities in units of 
a 0 (atomic units)- The results for the second dielectric virial 
coefficient, B £ , and the second Kerr virial coefficient, , 
however, are reported In units of cm 6 /mole 2 , and electrostatic 
units (esu) respectively (see Table V.l). The conversion factors 
which are therefore required are obtained in this appendix. 

To convert the units of B £ , Eq. (1.39) is quite straightforward. 
We have. 


B e = J dR R 2 a«(R) g(R) , 

where we assume that the integral is computed in atomic units (au), 
which, in this case, means that it has units of a® . Using 
N = 6.022169 x 10 23 /mole and la 0 = .529177 A = .529177 x 10" 8 cm, 
we have immediately that, 

B (cm 6 /mole 2 ) * .210 B (a®/mo1e 2 ) . 

£ £ ^ 

The conversion of to esu's is also straightforv/ard. To 
begin, we recall Eq. (1.45), 
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o t ,2m 2 f® 

B K " 4U5IT J 0 dR RZ 9( R ) » 

and note that, ir, general, the units of B K are length 9 /mole 2 -energy . 
One esu, then corresponds to taking the cm as the unit of length and 
one erg as the unit of energy. Taking k * K3S0622 x 10“ 16 erg/°K 
and noting that the atomic units of the integral are aj , it 
immediately follows that 6^ is obtained in esu's by multiplying 
the integral (in atomic units) by the factor: 


C = (16.61399/T) x 1Q“ 13 cmVaQ-erg-mole 2 . 


The particular values of C for T -• 4°K and 322°K are 
respectively: 


C(4°K) 


4=153497 x TO -13 


cm 9 

a§-erg-mole 2 




and 


C(322°K) 


>051 5962 x IQ" 13 


arr 

aS-erg-mole 2 


Throughout our work, the values of the various physical 
constants were taken from reference (83), 
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